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ABSTRACT 

This  work  is  an  examination  of  the  state-of-the-art  of  the 
two  layer,  integrated  type  atmospheric  mixed- Layer  models  with 
emphasis  on  applications  to  the  marine  regime.   A  linked  synoptic 
scale  mixed-layer  microphysical  scale  model  structure  is  suggested 
for  short  term  (24  hour)  forecasts  of  atmospheric  variables  with 
naval  applications.   Several  potential  deficiencies  in  the  present 
knowledge  are  identified  for  physical  parameterizations  (particu- 
larly the  entrainment  rate)  and  the  synoptic  scale  forcing 
requirements.   It  is  concluded  that  these  models,  though 
promising,  have  not  been  verified  for  a  wide  range  cf  atmospheric 
conditions . 


I  .    SUMMARY 

A.  INTRODUCTION 

It  is  clear  that  boundary  layer  models  hold  great  potential 
for  extending  and  expanding  local  atmospheric  assessments.   The 
integrated  type  mixed  layer  model  is  one  obvious  candidate  for 
this  application,  particularly  over  the  ocean.   Second  order 
closure  models  are  also  likely  candidates,  but  because  of  their 
much  greater  computer  requirements,  they  will  probably  be 
operational  somewhat  further  in  the  future.   The  integrated  model 
employs  rate  equations  for  the  mean  quantities  only,  assumes 
linear  vertical  flux  profiles  and  obtains  entrainment  rates  and 
turbulence  by  empirical  methods.   The  second  order  closure  models 
are  considerably  more  complex.   Since  our  own  area  of  expertise  is 
with  the  integrated  model,  we  shall  not  discuss  second  order 
closure  further. 

The  central  text  of  this  paper  is  an  exploration  of  the 
state-of-the-art  for  mixed  layer  models  with  a  focus  on  the 
sources  of  uncertainty  and  areas  needing  development.   We  will 
also  discuss,  somewhat  more  briefly,  the  present  situation  in 
obtaining  operational  output  and  problems  associated  with  actually 
implementing  a  field  version  of  this  model. 

B.  CONCLUSIONS 

Barring  some  very  extensive  breakthroughs  in  satellite  remote 
sensing  technology,  local  measurements  (such  as  radiosondes)  and 
meteorological  observations  will  be  required  to  provide  a  spectrum 
of  information  of  atmospheric  effects  on  Naval  operations. 
Extending  these  measurements  through  6-to-24  hour  forecasts  will 
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probably  require  on-board  computer  based  atmospheric  models 
supported  by  large  scale  numerical  field  models  with  regional 
window  capability  (for  example,  the  US  Navy  NORAPS  system).   The 
two  or  three  layer  integrated  type  marine  mixed  layer  model  is  a 
very  good  choice  for  these  local  assessments  because  it  is  ideally 
suited  to  the  good  horizontal  homogeneity  and  near  dynamic 
equilibrium  situations  that  often  characterize  the  open  ocean. 
Thus,  the  synoptic  scale  model,  the  local  data,  the  integrated 
boundary  layer  model  and  the  operational  output  parameterization 
could  be  linked  together  to  produce  local  boundary  layer 
forecasts.   A  sample  flow  diagram  for  such  a  system  is  given  in 
Fig.  1. 

The  feasibility  of  a  linked  synoptic-boundary 
layer-micrometeorology  model  is  primarily  dependent  upon  three 
issues:   1)   can  a  synoptic  scale  model  (with  regional  resolution) 
provide  sufficiently  accurate  forecasts  of  the  synoptic  forcing, 
2)   given  accurate  synoptic  forcing,  can  the  integrated  model 
provide  sufficiently  accurate  forecasts  of  the  boundary  layer 
evolutions  and  3)   can  the  relevant  operation  information  be 
extracted  from  the  mixed  layer  model  structure?   This  paper  will 
concentrate  on  question  two,  with  a  brief  background  on  the  third 
question.   The  first  question  is  beyond  the  area  of  expertise  of 
the  authors  of  this  paper  although  we  have  examined  certain 
alternatives  to  the  synoptic  scale  model  (namely,  single  station 
assessment  techniques) . 

The  following  conclusions  will  be  discussed  in  the 
examination  of  the  mixed  layer  model: 


1)  The  basic  model  structure  is  obtained  from  simple 
boundary  layer  averages  of  conservative  (or  near  conservative) 
quantities.   The  evolution  of  the  boundary  layer  is  then  predicted 
from  rate  equations  derived  from  the  conservative  expression. 
There  are  a  few  special  cases  to  consider  (temperature  and  mean 
wind),  but  in  principle  the  model  structure  is  extremely  simple. 

2)  The  internal  structure  of  the  model  requires  the  parame- 
terization of  certain  boundary  layer  energy  transfer  processes. 
These  parameterizations  are  in  fairly  good  shape  for  surface 
fluxes,  cloud  microphysics  and  radiative  flux  interactions.   There 
is  still  some  work  to  be  done  in  these  areas,  but  most  of  the 
questions  will  probably  be  resolved  within  the  next  few  years. 
Conversely,  the  entrainment  parameterizations  are  not  adequate  to 
handle  several  well  known  entrainment  mechanisms  (inversion  wind 
shear,  cloud  top  instability  or  atmospheric  wave  phenomena,  for 
example).   Since  inversion  wind  shear  conditions  are  common,  this 
particular  process  needs  further  investigation. 

3)  The  synoptic  and  mesoscale  effects  (subsidence  and 
advection)  are  often  as  important  as  entrainment.   Synoptic  scale 
numerical  models  are  very  important  here  because  not  only  can  they 
provide  the  relevant  data,  but  also  they  are  sources  of  forecasts 
of  synoptic  scale  forcing.   Whether  the  present  (or  near  future) 
models  will  be  adequate  has  not  been  determined.   The  mesoscale 
gradient  method  (using  satellite  remote  sensed  data  perhaps)  and 
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single  station  methods  may  be  able  to  provide  "nowcasts"  of 
certain  synoptic  terms.   For  example,  the  sea-surface  temperature 
field  will  almost  certainly  be  required  to  use  the  model 
successfully.   The  single  station  rawindsonde  can  provide  local 
values  of  subsidence  and  thermal  advection  but  the  accuracy  is 
marginal.   In  short,  the  required  large  scale  forcing  for  a 
complete  model  application  is  questionable. 

4)  The  algorithms  required  to  extract  operational  variables 
from  mixed  layer  parameters  are  in  an  advanced  state  of 
development.   There  are  some  problems  remaining  (continental 
aerosols,  for  example)  but  it  is  our  judgement  that  most  of  the 
bulk  atmospheric  characterizations  are  now  available  to  reasonable 
accuracy. 

5)  The  full  model  structure  as  defined  in  Fig.  1  is 
completely  unverified.   The  boundary  layer  model  has  only  been 
successfully  verified  on  a  few  selected  sets  of  atmospheric  data. 
Of  course,  most  of  the  individual  components  of  the  model  can  be 
verified  separately.   For  example,  a  candidate  radiation  flux 
profile  parameterization  can  be  tested  against  measured  radiation 
profiles  without  being  concerned  with  the  evolution  of  the 
boundary  layer.   Assuming  that  reasonably  accurate  boundary  layer 
expressions  and  large  scale  forcing  data  are  used,  then  testing 
the  model  against  data  is  basically  a  verification  of  the  model 
simplifying  assumptions  and  the  particular  entrainment  parameteri- 
zation chosen  by  the  researcher.   A  comprehensive  study  of  the 
accuracy  of  the  boundary  layer  model  over  a  wide  range  of  marine 
atmospheric  ccnditions  has  yet  to  be  performed. 
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II.   MODEL  BACKGROUND 

Boundary  layer  models  consider  two  basic  atmospheric  regions 
the  boundary  layer  ( Z <h  where  Z  is  the  height  above  the  surface) 
and  the  free  atmosphere  (Z>h)  separated  by  a  thin  transition 
region  at  height,  h.   The  boundary  layer  is  turbulent  while  the 
free  atmosphere  is  relatively  non-turbulent.   The  boundary  layer 
can  be  thought  of  as  an  interaction  region  between  two  large 
masses  -  the  ocean  and  the  free  atmosphere.   The  boundary  layer 
interacts  with  the  ocean  at  the  sea  surface  and  with  the  free 
atmosphere  at  the  transition  zone  (usually  referred  to  as  the 
inversion).   Because  of  the  strong  turbulence  in  the  boundary 
layer,  it  is  quite  homogeneously  mixed  vertically  (thus  the 
terminology  "well-mixed"  and  "mixed-layer"). 

The  well— mixed  nature  of  the  boundary  layer  implies  that 
certain  meteorological  variables  are  height  independent  for  Z <h . 
For  example,  the  water  vapor  mixing  ratio,  q  =  Q/p,  where  Q  is  th 
density  of  atmospheric  water  vapor  and  p  is  the  density  of  air 
(Fig.  2).   Thus  ,  we  assert  that  any  well-mixed  meteorological 
varable,  $  ,  can  be  represented  as 

4>  (Z)  =  $m  Z<h  (la) 

4>  (Z)  =  4>m  +  A$  +T  (z-h)  Z>h  (lb) 

where  tj>m  is  the  well-mixed  value,  &$   is  the  "jump"  at  the 
inversion  and  r  is  the  mean  vertical  gradient  of  t>   (assumed  zero 
for  Z<h)  in  the  free  atmosphere.   The  height  dependent  variable, 
4>(Z),  is  now  reduced  to  three  parameters:  p   m,  A<j>   and  T   .   This 
is  known  as  a  "zero-order"  model  since  we  allow  a  discontinuity  a 
the  inversion.   The  purpose  of  the  model  is  simply  to  predict  the 
time  evolution  of  these  three  parameters  for  each  variable  of 
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Figure  2.   Schematic  representation  of  the  two  layer  model 
structure  for  clear  sky  conserved  atmospheric 
variables-water  vapor  mixing  ratio,  Q/P,  and  virtual 
potential  temperature,  9y.   The  height  of  the 
well-mixed  layer  is  h   capped  by  an  inversion  of 
thickness  c*h . 


interest  and  relate  these  parameters  to  quantities  of  direct 
operational  relevance.   The  basic  method  is  to  relate  9  at  some 
time,  t,  {  Ht) )    to  an  initial  value,  ? i,  through  the  rate 
equation,  d$/dt,  for  a  small  time  step  t  -  tj_, 

d>(t)  =  4>i  +  (d$/dt)  (t  -  ti)  (2) 

The  rate  equations  are  relationships  of  d<f>/dt  to  the  physical 
processes  in  the  atmosphere.   For  the  boundary  layer,  the  major 
processes  are  interactions  with  the  surface  and  the  free 
atmosphere. 

Interaction  of  the  boundary  layer  air  with  the  ocean  surface 
is  the  primary  source  of  boundary  layer  turbulence.   (Note, 
however,  that  cloud  radiative  cooling  and  inversion  wind  shear  ca 
also  generate  turbulence.)   The  free  atmosphere  is  effectively 
"cut-off"  from  this  turbulence  by  the  inversion  zone.    The 
interaction  with  the  surface  is  characterized  by  surface  fluxes  o 
momentum  (wind),  sensible  heat  (temperature)  and  latent  heat  of 
evaporation  (water  vapor).   Since  the  fluxes  are  essentially  zero 
for  Z>h,  they  must  vary  with  height  within  the  boundary  layer.   I 
other  words,  the  vertical  flux  of  the  quantity  cp  ,  designated 
Fx(Z),  is  function  of  Z.   The  height  dependence  of  F   is  a 
critical  factor  in  the  temporal  evolution  of  1>  (Z).   For  a 
conservative  quantity,  the  rate  of  change  of  9  with  time  due  to 
the  flux  is  equal  to  the  vertical  gradient  of  the  flux 

d<t>/dt  =  -  dF^  (Z)/dZ  (3) 

The  well-mixed  assumption  (d$/dZ  =  0)  in  the  boundary  layer  is 
therefore  equivalent  to  assuming  that  dF^  (Z)/dZ  is  constant,  or 
that  Fx(Z)  has  a  linear  height  dependence. 


Equation  3  applies  to  changes  in  the  atmospheric  reference 
frame  ( Lagrangian) .   If  we  consider  D  at  a  point  fixed  above  the 
earth's  surface  (rather  than  a  point  that  translates  with  the 
atmospheric  motion) ,  then  we  must  account  for  the  possibility  that 
<J>  may  change  with  time  due  to  horizontal  gradients  in  <j>.   If  we 
define  the  mean  horizontal  wind  vector,  U,  the  horizontal 
gradient,  ^H  ,  the  mean  vertical  air  velocity,  W  and  the 
vertical  gradient  of  0,  d^/dZ,  then 

d<J>/dt-»-d<J>/dt   +   U-VH4>    +  "w   d^/dZ  (4) 

The  horizontal  term  is  referred  to  as  "advection"  and  the  vertical 
term  is  referred  to  as  "subsidence".   Thus,  the  equation 

d<J>/dt   =    -U-VH<j>    -  W    d<t>/dZ    -    dF^/dZ  (5) 

allows  us,  in  principal,  to  predict  the  evolution  of  a 

conservative  quantity  from  an  initial  measurement.   Note  that  most 

of  the  terms  in  Equation  5  are  time  and  height  dependent. 

The  next  step  in  applying  Equation  5  is  to  consider  the 

boundary  layer  and  free  atmosphere  separately  (in  the  spirit  of 

Equation  1).   In  the  boundary  layer,  we  integrate  Equation  5  to 

form  a  layer  average  (thus  the  terminology  "integrated"  model), 
h 

<X>  =  1/h  /  X(Z)dZ, 
o 


1/h  o/hdcj>m/dt   dZ    =    -1/h    /h[U-VH(t)    +   dF    /dZ]dZ  (6 


where   d?/dZ    =    0    and  4>    =  4>m   for    Z<h.      Therefore, 


d^/dt  =  (F^Q  -  F^h)/h  -  1/hJ  -U-VH*dZ  (7) 

o 

where  F^0  is  the  surface  flux  and  F,^  is  the  flux  at  the 
inversion. 

The  interaction  of  the  boundary  layer  with  the  surface  causes 
t>m  to  approach  the  characteristics  of  the  sea  surface  (via  the 
surface  flux,  F*0).   Similarly,  the  boundary  layer  interacts 
with  the  free  atmosphere  at  its  upper  boundary  through  a  process 
called  "entrainment "  (via  the  flux  F,^).   Entrainment  is  the 
process  whereby  free  atmospheric  air  is  brought  into  the  boundary 
layer  by  the  turbulent  erosion  of  the  transition  region.   Thus, 
entrainment  tends  to  increase  h  at  a  rate  characterized  by  the 
entrainment  velocity,  Wg,  and  causes  the  boundary  layer  to 
approach  the  properties  of  the  free  atmosphere.   The  increase  in  h 
is  expressed  by 

dh/dt  =  We  +  W  -  U  •  VRh  (8) 

Entrainment   produces    the    flux    at    the    top  of    the   boundary    layer. 
This    flux    is    analagous    to    the    surface    flux    at    the    sea-air    boundary 
(only    in   this    case    the    "surface"    is    the    transition    zone    to    the 
free    atmosphere) .       The    entrainment    velocity    is    simply    related    to 
F<t>h   by 


F<ph   "    "We  A*  O: 
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The    rate    equations    for    conservative    quantities    can    now   be 
written    for    both    atmospheric    regions 


d<j)m/dt    -    (F^Q    +   WeA±>)/h    -    <U-7H4>>  ;       z<h 


10a) 


d<|)/dt  =  -U  •  7H<(>  -  W  r.  ;   Z>h  (10b) 

d(A(J>)/dt  =  d<J)/dt)z=h+  -d<J>m/dt  -  I^dh/dt    ;   Z=h  (10c) 

These  equations  describe  only  the  evolution  due  to  turbulent 
transport  and  large  scale  advection.   For  atmospheric  quantities 
of  interest  (for  example,  temperature),  other  processes  can  appear 
as  source  and  sink  terms.   In  the  case  of  water  vapor,  vapor  can 
be  lost  due  to  formation  of  liquid  water  droplets  (condensation) 
and  subsequent  loss  from  rainout .   For  temperature,  there  is 
adiabatic  cooling  with  increasing  altitude,  latent  heat  released 
due  to  condensation  and  heat  transfers  associated  with  absorption 
and  emission  of  electromagnetic  radiation.   Such  effects  are  most 
easily  accounted  for  by  creating  an  appropriate  variable  which  is 
conserved.   For  example,  rather  than  apply  Equation  10  to  the 
water  vapor  mixing  ratio,  qv,  we  use  the  total  water 
qt  =  qv  +  q «,  where  q^  is  the  liquid  water  mixing  ratio.   For 
temperature,  we  define  the  moist  static  energy,  H, 

H  =  cp(T  +  .0098Z)  +  Leqt  (11) 

where  c_>  is  the  heat  capacity  of  air  and  Le  is  the  latent  heat 
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of  vaporization  of  liquid  water.   The  modifications  required  to 
account  for  radiative  flux  divergence  are  now  quite  simple;  rather 
than  use  FH ( Z )  in  Equation  10a,  we  use  the  sum  of  F^(Z)  and 
the  net  radiative  flux  R(Z). 

dHm/dt  =  (FHo  +  RQ  +  WeAH  -  RH)/h  -    <U- 7RH>  (12) 

The  application  of  the  total  model  now  requires  several 
steps:  1)  parameterization  of  physical  processes  (for  example, 
entrainment) ,  2)  development  of  techniques  to  obtain  the  relevant 
physical  data  (for  example,  subsidence  velocity  W  and  the 
horizontal  gradient  terms),  3)  development  of  the  relationships  of 
conserved  quantities  to  more  meaningful  meteorological  data  (for 
example,  the  relationship  of  temperature  to  moist  static  energy), 
and  4)  development  of  methods  to  calculate  operational  variables 
(for  example,  profiles  of  refractive-index  structure  function 
parameter,  Cn2 )  from  model  variables. 
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III.   MODEL  PERSPECTIVE 

Before  a  detailed  analysis  of  the  integrated  model,  it  is 
important  to  clearly  state  why  this  approach  cannot  be  a  cure-all 
for  all  atmospheric  predictive  problems.   This  model  contains 
certain  simplifying  assumptions  that  automatically  reduce  its 
generality  and  limit  the  scope  of  its  application.   Requirements 
for  the  two  layer  model  are: 

1)  Clearly  defined  mixed  layer  (or  at  least  a  capping 
inversion) . 

2 )  Surface  evaporation  leading  to  stratus  type  clouds 
confined  to  hc<Z<h  where  hc  is  the  lifting 
condensation  level. 

3)  Primary  source  of  vertical  transport  is  turbulent  mixing. 

4)  Moderate  or  predictable  synoptic  scale  thermal  changes  (a 
few  °C/day) . 

5)  Near  adiabatic  liquid  water  profiles  in  the  stratus 
cloud. 

6)  Moderate  clear-sky  net  radiative  heating  changes 
(+l°C/day). 

7 )  Moderate  or  predictable  horizontal  wind  divergence 
(|  D  |  ~  10~5  sec"1)  . 

Although  these  assumptions  are  quite  reasonable  over  a 
considerable  fraction  of  the  world's  oceans,  there  are  significant 
global  areas  where  they  break  down.   For  example,  in  tropical 
regions  much  of  the  vertical  transport  is  due  to  large  scale 
convergence  (the  Intertropic  Convergence  Zone).   Because  the  total 
water  vapor  content  of  the  air  is  greater  (higher  sea  surface 
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temperatures)  the  clouds  have  much  more  buoyancy  from  latent  heat 
of  condensation.   These  highly  buoyant  clouds  cannot  be  confined 
by  the  inversion  so  they  break  through  and  rise  to  great  altitudes 
as  classic  free  convection  cells.   The  large  bouyant  and  turbulent 
mixing  processes  in  the  cloud  tend  to  destroy  the  adiabatic 
(equilibrium)  liquid  water  profile.   Rather  than  entraining  upper 
layer  air  into  the  mixed  layer,  the  convective  cells  draw  air  out 
of  the  mixed  layer  into  the  cell  (and,  therefore,  the  upper  air) 
causing  a  lowering  of  the  mixed  layer  height.   This  process  is 
called  cloud  induced  subsidence  or  detrainment.   Considerable  work 
is  now  underway  to  apply  the  integrated  model  to  this  case  using 
an  additional  layer  (two  inversions)  and  parameterizations  of 
cumulus  clouds. 

Another  case  for  which  this  model  cannot  be  used  is  near 
frontal  zones.   Clearly,  almost  every  assumption  listed  can  be 
expected  to  fail  in  the  vicinity  of  a  strong  front.   As  a  result, 
we  do  not  expect  the  model  to  be  very  useful  at  higher  latitudes 
during  the  winter. 

Despite  these  disclaimers,  the  model  does  have  a  wide  range 
of  utility.   A  recent  paper  by  Kraus  and  Schaller  shows  examples 
of  mixed  layer  model  application  to  four  different  climate 
regions:  tradewind,  California  coast,  Norwegian  Sea  and  the  Arctic 
Ocean  (summer) . 
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IV.   PHYSICAL  PROCESSES 

Parameterization  of  the  physical  processes  is  essentially  the 
art  of  calculating  the  unknown  from  the  known.   For  mixed  -layer 
models,  the  known  is  primarily  the  mean  meteorological  profiles, 
basic  atmospheric  thermodynamics  and  a  compendium  of  relevant 
atmospheric  data  (such  as  the  radiant  intensity  from  the  sun  at 
the  top  of  the  atmosphere).   This  section  will  be  devoted  to  an 
examination  of  various  parameterizations  used  to  obtain  each  term 
in  the  rate  equations. 

A.   Surface  Fluxes 

The  bulk  aerodynamic  method  is  used  to  obtain  the  surface 
fluxes  of  momentum,  sensible  heat  and  water  vapor.   The  basic 
assumption  of  this  method  is  that  the  flux  is  driven  by  the 
sea-air  bulk  difference  of  the  variable 
V  -  cl/2  (<DS  -  >m)u*  (13) 

where  the  subscript  "s"  refers  to  the  sea-surface  value  and  U* 
is  the  friction  velocity 

U*  =  (Fuo/p)l/2  =  CU1/2(J  (14) 

where  U  is  the  wind  speed  near  the  surface  (Z^lOm).   The  quantity 
C,  is  known  as  the  drag  coefficient  (typically  C  1/2  =  0.035). 
The  bulk  method  has  been  so  well  verified  that  it  is 
non-controversial  and  further  discussion  in  the  mixed  layer  model 
context  is  not  needed. 
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B.   Clouds 

There  are  several  factors  to  consider  when  dealing  with 
boundary  layer  clouds:  cloud  thickness,  liquid  water  profile  and 
cloud  droplet  size  distribution. 

1 )   The  most  common  assumption  is  that  the  cloud  bottom  is  at 
a  height  hc  equal  to  the  lifting  condensation  level.   Since  we 
are  assuming  constant  potential  temperature 

(9=T  +  .0098Z)  and  water  vapor  mixing  ratio  (qv  =  qt)  in  the 
noncloud  region  (Z<hc),  hc  is  easily  calculated  from  the  near 
surface  temperature  (T)  and  dew— point  temperature  (Td) 

hc  =  (T  -  Td)/0. 00804  (15) 

where  0.00804  is  the  difference  in  the  lapse  rate  of  T  and  Td 
under  well  mixed  conditions.   However,  because  hc  is  very 
sensitive  to  the  9  and  q  profile  in  the  mixed  layer,  small 
deviations  from  the  well-mixed  assumption  (which  for  most  other 
purposes  are  not  significant)  cause  large  changes  in  hc  (and, 
therefore,  the  cloud  radiative  effects).   Since  warm,  dry  air  is 
entrained  into  the  boundary  layer  at  the  top,  it  is  generally 
observed  that  qtm  decreases  and  9m  increases  slightly  with 
increasing  height  in  the  mixed  layer.   Therefore,  hc  is  actually 
considerably  greater  (thinner  cloud)  than  hc  calculated  from 
near  surface  observations  using  Equation  15.   The  error  can  often 
be  more  than  100m,  which  is  significant  for  clouds  that  are 
typically  200  to  500m  thick.   At  the  moment,  only  "ad  hoc"  methods 
are  available  for  correcting  for  this  effect.   These  methods  allow 
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qtm  and  Hm  to  have  small  gradients  that  depend  on  the 
surface  flux  and  the  entrainment  rate. 

2)  Most  mixed  layer  models  calculate  the  liquid  water 
profile  in  the  cloud  based  on  the  same  adiabatic  assumption 
(qtm  and  Hm  are  height  independent  for  Z <h ) .   In  other 

words,  given  qt,  H  and  Z  one  can  uniquely  determine  q, . 
Experimentally,  it  has  been  found  that  the  liquid  water  profile  is 
almost  always  less  than  the  adiabatic  value  for  the  reasons 
explained  in  part  1.   If  we  use  a  model  that  allows  a  slight 
decrease  of  qtm  and  a  slight  increase  of  Hm,  then  the  cloud 
liquid  water  profile  will  be  decreased.   Since  the  model  contains 
arbitrary  constants,  we  simply  make  adjustments  to  obtain 
agreement  with  the  measurements  (roughly  a  30%  reduction  of  liquid 
water  from  the  adiabatic  value). 

3 )  Cloud  droplet  spectra  used  in  most  mixed  layer  models  are 
based  on  strictly  empirical  parameterizations  of  field 
measurements.   Unfortunately,  the  measurements  are  extremely 
sparse  and  not  very  consistent.   The  most  credible  data  indicates 
a  fairly  constant  total  number  density,  N,  with  the  mean  radius  of 
the  cloud  drops  increasing  with  height  in  the  cloud  (rm  —  2ym  to 
rm  —  15]jm)  .   Sample  data  are  shown  in  Fig.  3.   There  are  data 
that  indicate  fairly  constant  rm  but  an  increasing  N  with 
altitude  (for  example,  100  cm~3  to  300  cm"^).   There  is 
considerable  evidence  that  near  the  cloud  top,  variations  in 
liquid  water  are  due  to  variations  in  N  while  rm  remains  fairly 
constant.   It  has  been  suggested  that  this  is  due  to  the 
inhomogeneous  mixing  of  cloud  air  with  newly  entrained  clear  air. 
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Figure  3.   Cloud  droplet  spectral  parameters  from  mid-Atlantic 

stratus  (Ref.13).   The  dispersion  is  the  ratio  of  the 
standard  deviation  to  the  mean  radius  of  the  size 
spectrum. 


Obviously,  the  integrated  cloud  droplet  volume  spectrum  should 
give  the  correct  liquid  water  at  each  altitude 

oo 

q^p  =  Q  (Z)  =  pw  f  4/3  Tir3(dN/dr)dr  (16) 

where  pw  is  the  density  of  water,  r  the  droplet  radius  and  dN/dr 
the  droplet  size  spectrum.   There  is  also  a  scale  problem  because 
most  of  the  liquid  water  is  in  the  larger  drops  (r>10pm)  while 
most  of  the  shortwave  (solar)  radiation  interaction  is  determined 
by  smaller  droplets.   Therefore,  one  can  have  two  different 
empirical  droplet  distributions,  both  scaled  by  the  liquid  water, 
with  quite  different  short  wave  radiation  absorption.   In  fact,  it 
has  been  shown  that  the  longwave  emission,  short  wave  absorption 
and  the  entrainment  rate  are  all  dependent  upon  the  droplet 
distribution.  • Since  this  is  a  field  of  very  active  research,  it 
is  probable  that  good  empirical  models  will  be  available  very 
soon. 

C.   Radiation 

Mathematical  modeling  of  radiative  flux  transfer  is  an 
extensively  developed  subject.   In  general,  the  radiative  effects 
can  be  calculated  to  a  precision  far  exceeding  the  quality  of  the 
input  data.   In  fact,  the  main  problem  is  to  choose  a  method  of 
minimum  complexity  appropriate  to  the  accuracy  of  the  data  to  be 
used.   Unfortunately,  even  the  simple  radiative  transfer  models 
are  still  extremely  complex  in  comparison  to  most  other  physical 
parameteri-ations  used  in  the  mixed  layer  model.   Uncertainty  in 
background  aerosols,  atmospheric  absorbing  gases  (water  vapor, 
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carbon  dioxide,  ozone)  and  cloud  droplet  size  spectra  are  obvious 
sources  of  error  in  radiation  calculations.   A  comparison  of  two 
models  with  aircraft  measurements  of  the  radiative  flux  divergence 
is  shown  in  Fig.  4. 

Most  models  consider  longwave  and  shortwave  radiation 
separately.   Typically,  radiation  effects  are  calculated  only 
inside  the  cloud  layer  (hc  <  Z<h).   Radiative  heating  and 
cooling  rates  in  clear  air  are  on  the  order  of  l°C/day  which  is 
approximately  the  accuracy  of  even  the  best  models,  given  typical 
input  data.   Earlier  models  that  simply  assign  a  constant  long 
wave  cooling  rate  to  the  cloud  are  now  considered  to  be  too  crude 
because  they  neglect  the  dependence  of  the  downward  flux 
(effective  sky  temperature)  on  upper  atmosphere  properties 
(primarily  water  vapor  content)  and  they  assume  the  cloud  is 
totally  "black"  regardless  of  cloud  thickness.   Actually,  the 
emissivity,  E,  of  clouds  in  the  long  wave  region  does  not  approach 
unity  until  h-hc  ~  100m.   In  the  longwave  region,  the  absorption 
is  not  extremely  sensitive  to  variations  in  the  droplet  spectra 
typically  found  in  clouds;  it  turns  out  that  the  emissivity  can  be 
parameterized  in  terms  of  the  total  liquid  water,  W 


W  =     pq£dZ  (17a) 


E  =  1  -  exp(-0.138W)  (17b) 
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Obviously,  the  integrated  cloud  droplet  volume  spectrum  should 
give  the  correct  liquid  water  at  each  altitude 

q^p  =  Q  (Z)  =  pw  j     4/3  7rr3(dN/dr)dr  (16) 

where  pw  is  the  density  of  water,  r  the  droplet  radius  and  dN/dr 
the  droplet  size  spectrum.   There  is  also  a  scale  problem  because 
most  of  the  liquid  water  is  in  the  larger  drops  (r>10um)  while 
most  of  the  shortwave  (solar)  radiation  interaction  is  determined 
by  smaller  droplets.   Therefore,  one  can  have  two  different 
empirical  droplet  distributions,  both  scaled  by  the  liquid  water, 
with  quite  different  short  wave  radiation  absorption.   In  fact,  it 
has  been  shown  that  the  longwave  emission,  short  wave  absorption 
and  the  entrainment  rate  are  all  dependent  upon  the  droplet 
distribution.   Since  this  is  a  field  of  very  active  research,  it 
is  probable  that  good  empirical  models  will  be  available  very 
soon. 

C.   Radiation 

Mathematical  modeling  of  radiative  flux  transfer  is  an 
extensively  developed  subject.   In  general,  the  radiative  effects 
can  be  calculated  to  a  precision  far  exceeding  the  quality  of  the 
input  data.   In  fact,  the  main  problem  is  to  choose  a  method  of 
minimum  complexity  appropriate  to  the  accuracy  of  the  data  to  be 
used.   Unfortunately,  even  the  simple  radiative  transfer  models 
are  still  extremely  complex  in  comparison  to  most  other  physical 
parameterizations  used  in  the  mixed  layer  model.   Uncertainty  in 
background  aerosols,  atmospheric  absorbing  gases  (water  vapor, 
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carbon  dioxide,  ozone)  and  cloud  droplet  size  spectra  are  obvious 
sources  of  error  in  radiation  calculations.   A  comparison  of  two 
models  with  aircraft  measurements  of  the  radiative  flux  divergence 
is  shown  in  Fig.  4. 

Most  models  consider  longwave  and  shortwave  radiation 
separately.   Typically,  radiation  effects  are  calculated  only 
inside  the  cloud  layer  (hc  <  Z<h).   Radiative  heating  and 
cooling  rates  in  clear  air  are  on  the  order  of  l°C/day  which  is 
approximately  the  accuracy  of  even  the  best  models,  given  typical 
input  data.   Earlier  models  that  simply  assign  a  constant  long 
wave  cooling  rate  to  the  cloud  are  now  considered  to  be  too  crude 
because  they  neglect  the  dependence  of  the  downward  flux 
(effective  sky  temperature)  on  upper  atmosphere  properties 
(primarily  water  vapor  content)  and  they  assume  the  cloud  is 
totally  "black"  regardless  of  cloud  thickness.   Actually,  the 
emissivity,  E,  of  clouds  in  the  long  wave  region  does  not  approach 
unity  until  h-hc  -  100m.   In  the  longwave  region,  the  absorption 
is  not  extremely  sensitive  to  variations  in  the  droplet  spectra 
typically  found  in  clouds;  it  turns  out  that  the  emissivity  can  be 
parameterized  in  terms  of  the  total  liquid  water,  W 

h 
W  =  f   pq£dZ  (17a) 

c 

E  =  1  -  exp(-0.138W)  (17b) 
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Figure  4.   Net  shortwave  (darkened  boxes)  and  longwave  (clean 

boxes)  heating  rates  (flux  divergences)  measured  from 
an  aircraft  are  shown  on  the  left.   The  right  side 
depicts  the  total  radiation  heating  rate  profile.   The 
insert  compares  these  measurements  with  two  recent 
models  (Ref.13). 


Figure  5.   The  entrainment  rate,  We,  as  a  function  of  total 

integrated  cloud  liquid  water  content,  W  (Eq.  16a). 
The  influence  of  solar  radiation  is  shown  for  three 
different  solar  zenith  angles  (Ref.  7).   Note  8  =  90 
correspond  to  zero  solar  flux. 
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The    cloud    top    cooling    rate    can   be    estimated    from    the    cloud    top 
temperature,    T^ ,    and    the    effective    sky    temperature,    Tsky, 
using    the    Stef an-Bolzmann    law 

ARh    =    Ea(Th4-Tsky4) 

where  a  is  Stefan's  constant.   Similarly,  the  cloud  bottom  is 
heated  by  infrared  radiation  from  the  sea  surface  (which  is 
usually  warmer  than  the  cloud).   Although  this  approach  is  quite 
simple  for  calculating  the  cooling  rate,  it  leads  to  complications 
when  calculating  the  effect  the  radiative  flux  divergence  has  on 
the  entrainment  rate  (which  is  sensitive  to  the  actual  vertical 
distribution  of  flux  divergence  within  the  cloud).   This  will  be 
discussed  further  in  the  section  on  entrainment  rate. 

In  the  shortwave  region,  the  emissivity  approach  (Equation 
16)  is  not  very  useful  because  shortwave  extinction  is  dominated 
by  scattering  (as  opposed  to  absorption).   Thus,  the  direct  solar 
component  incident  at  the  cloud  top  is  attenuated  rapidly  as  the 
cloud  is  penetrated,  but  the  vast  majority  of  the  radiation  lost 
from  the  direct  beam  does  not  warm  the  cloud,  but  is  simply 
scattered.   This  scattered  radiation  forms  a  second  component 
usually  referred  to  as  the  "diffuse"  solar  radiation.   A  small 
fraction  of  the  direct  and  diffuse  radiation  is  absorbed  by  the 
cloud  droplets,  thus  warming  the  boundary  layer.   The  amount  of 
heating  depends  on  the  solar  zenith  angle,  the  droplet  spectrum, 
the  droplet  index  of  refraction  and  the  albedo  of  the  sea  surface. 
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A  number  of  methods  are  available  for  calculating  the  short 
wave  heating  in  clouds.   The  so-called  "two  stream"  approximation, 
where  the  scattered  radiation  is  divided  into  a  forward  component 
and  a  component  with  an  angular  distribution  calculated  from  the 
droplet  spectrum,  is  certainly  accurate  enough  for  mixed  layer 
model  applications.   This  method  divides  the  cloud  into  several 
(typically  four)  layers  in  the  vertical  to  properly  account  for 
the  redistribution  of  direct  and  diffuse  radiation  as  a  function 
of  height  in  the  cloud.   An  excellent  review  of  this  subject  has 
been  published  by  Zdunkowski  et  al  (1980). 
D.   Entrainment 

Fluxes  at  the  top  of  the  boundary  layer  are  due  to  the 
entrainment  of  free  atmospheric  air.   We  must  determine  this  flux 
for  each  variable  of  interest;  for  N  variables  we  have  N  unknown 
fluxes,  F,^.   Eq .  9  relates  F  h  to  the  entrainment  velocity, 
We,  and  the  jump  in  <£  across  the  inversion,  &$.       It  appears  we 
have  traded  one  unknown  (F  ^)  for  two  (We  and  Ad>).   Actually, 
one  value  of  We  is  appropriate  for  all  F  h  and  the  jumps  are 
not  unknowns  which  must  be  parameterized  but  instead  are 
properties  which  can  be  modeled  from  the  initial  measurements. 
Thus,  we  have  reduced  the  N  unknown  fluxes  to  one  unknown 
parameter,  WQ .   The  problem  left  is  to  determine  We . 

The  entrainment  process  has  been  investigated  through  labora- 
tory simulations,  three-dimensional  numerical  models,  field  meas- 
urements (primarily  with  aircraft)  and  theoretical  simplification 
of  the  turbulent  kinetic  energy  balance  expressions  (this  is  the 
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so-called  "closure"  assumption) .   The  number  of  entrainment 
relations  available  is  approxiamtely  equal  to  the  number  of 
investigators  in  the  field,  which  is  an  indication  of  the 
complexity  of  the  problem.   Basically,  since  entrainment  is  the 
erosion  of  the  inversion  layer  by  boundary  layer  turbulence,  the 
entrainment  rate  depends  on  the  intensity  of  turbulence  at  the  to 
of  the  boundary  layer,  the  stability  of  the  inversion  layer  and 
the  stability  of  the  upper  atmosphere.   The  virtual  potential 
temperature,  9y,  is  an  index  of  the  hydrostatic  stability  of 
air.   If  one  parcel  of  air  has  greater  $v  than  another,  it  is 
said  to  be  more  buoyant.   Thus,  the  jump  in  ^v  across  the 
inversion,  ^9V,  is  a  primary  indicator  of  the  resistance  of  the 
inversion  to  turbulence.   Turbulence  at  the  inversion  available 
for  entrainment  must  be  either  produced  at  the  inversion  or 
transported  to  the  inversion  from  production  regions  in  the  lower 
parts  of  the  boundary  layer.   Thus,  all  processes  that  create 
turbulence  in  the  boundary  layer  are  potential  sources  of 
entrainment  energy.   Since  turbulence  is  rapidly  dissipated  by 
viscous  forces,  sources  near  the  inversion  tend  to  be  more 
efficient  at  enhancing  the  entrainment  process. 

The  simplest  case  to  deal  with  is  the  cloud  free  boundary 
layer  with  negligible  jump  in  the  mean  wind  vector  at  the 
inversion.   Since  the  boundary  layer  turbulence  in  this  case  is 
driven  by  the  surface  flux  only,  the  entrainment  flux  is  assumed 
to  be  a  fraction  of  the  surface  flux 

F9  h-  -WeA8v  =  -A  Fa  Q  (18) 


where  the  fraction,  A,  has  been  found  to  be  about  0.2.   Since 

A9V  must  be  positive  (otherwise,  no  inversion  exists),  the 

negative  value  of  F^  ^  is  simply  a  reflection  of  the  fact  that 

v 

the  boundary  layer  turbulence  must  do  work  to  entrain  free 
atmospheric  air  across  the  highly  stable  inversion.   Thus,  the 
buoyancy  entrainment  flux  represents  a  turbulence  energy  loss 
mechanism  which  (along  with  viscous  dissipation)  is  balanced  by 
turbulence  produced  in  the  lower  boundary  layer  that  has  been 
transported  to  the  inversion. 

Although  Eq .  18  is  elegant  in  its  simplicity,  it  is  only 
applicable  in  cloud  free,  light  wind  situations.   Clouds  are 
important  because  they  are  a  source  of  longwave  cooling  at  the  top 
of  the  boundary  layer.   Since  cooling  at  the  top  is  equivalent  to 
heating  at  the  bottom,  clouds  are  source  of  convective  (buoyant) 
production  of  turbulence  and,  therefore,  entrainment.   Since  the 
longwave  and  shortwave  fluxes  are  distributed  differently  in  the 
cloud,  the  entrainment  rate  also  depends  on  the  shortwave  profile 
(Fig.  5).   Assuming  one  can  calculate  the  cloud  properties  and 
and  radiative  flux  profiles,  then,  in  principle,  the  cloud 
entrainment  contribution  can  be  determined.   Some  controversy 
still  exists  on  the  correct  manner  to  relate  the  vertical 
structure  to  the  horizontal  averaging  method.   The  horizontal 
averaging  process  enters  the  problem  because  of  the  "bumpy"  nature 
of  the  top  of  the  boundary  layer.   It  is  not  known  what  fraction 
of  cloud  top  radiative  cooling  is  confined  to  the  mixed  layer. 
This  is  presently  an  active  research  area. 
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The  inversion  wind  shear,  which  occurs  if  there  is  a  jump  in 
the  mean  wind  velocity  (  AU),  has  not  received  much  attention. 
The  importance  of  wind  shear  for  a  given  situation  can  be 
estimated  by  the  following  expression 

We  =  Weo/(l-Ri-l)  (19) 

where  Weo  is  the  entrainment  rate  due  to  bouyancy  effects  and 
Rj_  is  the  inversion  Richardson  number 

R±  =  gah  A9v/(t|AU  |  2)  (20) 

where  g  is  the  gravitational  acceleration  and  ah  is  the  thickness 
of  the  entrainment  region.   If  Rj_>>1,  then  one  can  neglect 
inversion  wind  shear  effects.   If  R^  is  on  the  order  of  one, 
then  it  is  likely  that  conventional  entrainment  parameterizations 
will  not  be  adequate. 
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V.   Synoptic  Scale  Forcing 

The  role  of  synoptic  scale  air  motion  on  fixed  point  evolu- 
tions was  expressed  in  Eq.  5.   For  example,  typical  atmospheric 
horizontal  temperature  gradients  are  on  the  order  of  1°C  for  100km. 
At  |  U|  =  10  m/s,  this  could  result  in  an  advection  temperature 
change  as  large  as  3.6°C  in  12  hours.   From  the  model  application 
point  of  view,  one  regards  U  and  d<f/dZ  as  part  of  the  model 
initialization  and  rate  equation  structure  (Eq.  1  and  Eq.  10), 
just  as  the  A<£  data  were  regarded  in  the  entrainment  section.   The 
quantities  W  and  7H<f>  are  not  predicted  by  the  model  but  are 
parameters  that  must  be  specified  in  order  to  make  valid 
predictions.   There  are  three  methods  that  are  potential  sources 
of  this  information:   1)   synoptic  scale  numerical  models,  2) 
multiple  station  gradient  measurements  and  3 )  single  station 
techniques . 

The  mean  vertical  velocity,  (subsidence)  and  the  horizontal 
wind  gradients  (divergence)  are  phrases  that  tend  to  be  used 
interchangeably  because  they  are  directly  related  through  mass 
conservation 

V-  (U+Wp£)=Vh'U+  dW/dZ  =  D  +  dW"/dZ  =  0  (21) 

Therefore,  the  mean  vertical  motion  at  height  Z  is  simply 
Z 


W(Z)=-|      VH-UdZ 


(22; 

o 

Often    the    divergence    is    assumed    to    be    independent   of  height   below 

5km  giving 

W   =    -DZ  (23) 
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A.  Synoptic  Scale  Numerical  Models 

Synoptic  scale  methods  can  be  used  to  calculate  W  .   The  glo- 
bal pressure,  temperature  and  geopotential  fields  theoretically 
contain  the  information  required  to  calcualte  the  horizontal  wind 
field,  u*,  the  subsidence,  and  the  horizontal  gradients  of  U  and  T. 
Subsidence  can  be  obtained  from  direct  divergence  calculation, 
several  different  vorticity  methods  and  thermodynamic  methods.   At 
the  moment,  the  global  field  estimates  over  the  ocean  are  based  so 
heavily  on  data  obtained  over  land  that  these  methods  are  consi- 
dered unreliable  for  subsidence  and  gradient  calculations.   When 
over  ocean  data  become  routinely  available  with  reasonable  density 
and  vertical  resolution,  synoptic  methods  will  be  valuable  for  the 
wind  field,  subsidence  and  temperature  advection  terms  above  the 
boundary  layer.   There  are  presently  no  plans  to  develop  this 
resource  for  water  vapor  and  aerosol  advection. 

B.  Mesoscale  gradients 

Suppose  we  wish  to  consider  the  evolution  of  the  boundary 
layer  over  a  12  hour  period.   For  a  typical  wind  speed  (10  m/s), 
we  are  concerned  with  an  area  roughly  400km  on  a  side.   This  is 
really  a  more  mesoscale  area  of  interest.   If  we  have  available  a 
network  of  atmospheric  measurements  of  all  relevant  phenomena  from 
the  surface  to  several  km  altitude,  distributed  over  the  area  of 
interest,  then  the  horizontal  and  vertical  gradients  and  the 
subsidence  can  be  calculated.   The  classic  method  of  obtaining  the 
divergence  is  to  integrate  the  horizontal  wind  field  components 
perpendicular  (U  )  to  a  closed  horizontal  path 


7  H  •  U  =  (1/A)  6   UDd2,  (24) 


2  3 


where  A  is  the  area  enclosed  by  the  horizontal  path.   A  minimum  of 
three  separate  locations  (a  triangle)  is  required. 

A  more  general  method,  which  can  be  applied  to  a  horizontal 
field  of  measurements  of  any  variable,  is  to  perform  a  linear 
regression  on  all  data  points  with  horizontal  coordinates  X  and  Y 
as  variables  and  slope  parameters  4>  x  and  iy,  Z  held  constant 

<J>  (x,y)    =   a   +  <t>  XX    +  *yY  (25) 

Since    d^/dx    =    <J>X   and   d<}>/dy   =    <jy,    then 

U  •  VH0  =  ux4>x  +  Uycpy  (26) 

describes  the  horizontal  advection  of  the  quantity  <$> . 

If  sufficient  data  are  available  (admittedly  a  big  "if"), 
the  mesoscale  field  method  appears  to  be  ideal.   However,  in 
addition  to  the  great  difficulty  involved  in  operating  a  marine 
data  network,  measurement  inaccuracies  and  atmospheric  variability 
impose  inherent  limits  on  the  accuracy  of  this  method.   For 
example,  a  2  m/s  inaccuracy  in  wind  speed  (typical  for 
conventional  rawindsondes )  over  a  200km  triangle  represents  an 
uncertainty  in  divergence  of  t  1  x  10~5  s"1.   This  could 
cause  an  uncertainty  in  the  height  of  the  boundary  layer  of 
approximately  a  factor  of  two  within  24  hours .   The  problem  of 
atmospheric  variability  can  be  solved  by  averaging  a  number  of 
near  simultaneous  (within  one  hour)  atmospheric  profiles  at  each 
location.   For  a  one  minute  average,  typical  wind  velocity 
variability  is  about  15%  of  the  mean  with  larger  variability  at 
low  wind  speeds.   Thus,  wind  velocity  uncertainty  due  to 
variability  is,  again,  on  the  order  of  2  m/s  at  typical  mean  wind 
speed  for  a  single  sounding.   It  appears  that  the  mesoscale  mean 
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field  gradient  methods  must  await  the  deployment  of  a  more 
accurate,  greater  horizontal  density  ocean  measurement  system. 
C.   Single  Station  Analysis 

Since  local  single  station  data  will  almost  certainly  be 
required  for  the  model  initialization,  a  method  that  allows 
estimation  of  the  synoptic  terms  from  these  data  is  ideal.   In  the 
case  of  subsidence,  the  existing  methods  require  examination  of 
successive  radiosonde  profiles.   One  obvious  method  (basically  a 
"bootstrap"  method)  is  to  simply  adjust  the  value  of  W  used  in  the 
previous  time  period  to  obtain  agreement  of  the  model.   This 
assumes  that  the  model  entrainment  rates  and  boundary  layer  height 
advective  terms  are  correct  and  the  only  remaining  source  of  error 
in  the  prediction  of  h  is  the  subsidence  velocity.   One  must  also 
assume  that  the  average  subsidence  over  the  previous  time  period 
is  applicable  to  the  current  time  period. 

A  second  technique  utilizes  the  heating  of  the  upper  atmos- 
phere (Z>h)  under  subsidence  conditions.   At  a  given  height,  the 

change  of  potential  temperature,  d„,  over  the  time  interval,  dt, 

9 

is  related  to  the  subsidence  velocity  by  (Eq.  10b) 

W  =  (u  •  7H9  -  d9/dt)/Te  (27) 

If  advection  is  neglible  and  the  divergence  is  independent  of 

height,  then  (using  Eq.  23) 

D  -  d9/( rQ  dtZ)  (28) 

Note  that  under  these  conditions  the  heating  rate  is  proportional 

to  height,  therefore  that  lapse  rate,  r~,  will  not  remain  constant 

in  time.   Thus,  the  change  in  lapse  rate  with  time  can  be  used  to 

estimate  divergence 
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D  =  £n(rQi/rQ2)/dt  (29) 

where  ?q ,  and  ?q .are  the  lapse  rates  at  times  t2  and  t]_  =  t2  -   t. 

If  accurate  estimates  of  temperature  advection  are  available, 

subsidence  can  be  accurately  obtained  using  the  thermodynamics 

method  ( Eq .  27)  with  a  time  series  of  single  station  temperature 

profiles.   One  technique  is  available  to  obtain  single  station 

temperature  advection,  the  thermal  wind  method.   This  method 

utilizes  the  known  change  of  the  mean  wind  vector  as  a  function  of 

height  due  to  local  gradients  of  the  temperature  field.   The  wind 

vector  change  is  described  by  the  thermal  wind  equation 

t        =  U]_  -  uQ  =  R/f  *  (k  x  V9  )  J,n(P0/P1)  (30) 

where  1  and  0  designate  different  altitudes,  P  is  the  pressure,  R 

the  gas  constant,  and  f  =  2fi  sin( latitude)  is  the  coriolis 

parameter.   If  we  designate 

7  9=  d9/dx  i  +  d9/dy  ? ,    then 
H 

de/dx  =  +f(Vly  -  Voy)/(R  in(P0/P!))  (31a) 

dG/dy  =  -f(Vlx  -  Vox)/(R  ln(P<D/P1    ))  (31b) 

If  we  further  designate  U  =  (U]_  +  UQ  )  /2 ,  then  the  advection 

heating  is 

-U  •  7H0   =  f(V2   +  V2   -  V2.  -  V2  )/(2R£n(P  /P,  )  )  (32) 

"         lx    oy    ly    ox        o   1 

This  method  does  not  work  in  the  boundary  layer  where  wind 
gradients  are  also  caused  by  the  surface  stress . 

In  principle,  Eq .  32  and  Eq .  27  can  be  used  as  single  station 

■fc 
assessments  of  W  and  -U •  7H9 .   However,  since  these  are 

difference  methods  (either  between  different  soundings  or 

different  levels  within  a  sounding),  they  are  subject  to  the  same 

inherent  uncertainties  due  to  measurement  inaccuracy  and  normal 
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atmospheric  variability  described  in  the  section  on  the  gradient 
method.   Fig.   6  shows  a  successful  comparison  of  thermal  wind  and 
gradient  method  estimations  of  temperature  advection. 
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Figure  6.   Thermal  advection  heating  rate,  -U  •  79,  as  a 

function  of  altitude.   The  data  are  obtained  from  four 
rawindsonde  stations  in  the  southern  California  coastal 
region.   The  solid  line  is  calculated  using  the  thermal 
wind  equation  (Eq.32)  and  the  dashed  line  is  from  the 
gradient  method  (Eq.26).   The  close  agreement  of  the 
two  methods  is  not  typical. 
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VI.   USER  OUTPUT 

The  basic  product  of  the  integrated  model  is  a  prediction  of 
the  time  evolution  of  the  mean  profiles  (Eq.  4)  and  the  height  of 
the  boundary  layer.   Also  available  are  the  various  quantities 
required  to  drive  the  model  -  surface  fluxes,  cloud  thickness, 
liquid  water  profiles,  long  and  short  wave  radiation  flux  profile; 
plus  external  specifications  of  advective  and  subsidence  effects. 
This  information,  which  is  the  direct  output  of  the  model,  has 
numerous  applications.   Certain  other  relevant  data,  such  as  the 
electromagnetic  refractive- index  profile,  are  easily  calculated 
from  the  model  variables  (in  this  case,  the  temperature  and  water 
vapor  profiles).   There  is  a  third  class  of  information 
(turbulence  properties,  for  example)  which  require  empirical 
models  for  calculation  from  mixed  layer  mean  profiles. 

Typical  Navy  applications  would  include  basic  meteorology, 
optical  propagation,  EM  (radio  and  radar)  propagation,  atmospheric 
dispersion,  remote  sensing  and  oceanographic  forcing.   In  this 
group,  the  EM  propagation  (IREPS),  meteorological  and  oceanogra- 
phic model  structure  is  very  well  developed.   The  optical  propaga- 
tion, atmospheric  dispersion  and  remote  sensing  applications  are 
not  as  completely  developed  (they  probably  lag  the  development  of 
the  mixed  layer  evolutionary  model  at  this  stage) .   As  an  example, 
the  state-of-the-art  for  optical  propagation  is  discussed  in  an 
appended  publication.   Fig.  7  illustrates  how  the  relevant  physic;, 
quantity  (in  this  case,  the  aerosol  profile)  must  be  cast  in  the 
mixed  layer  model  format  to  utilize  the  model  predictions. 
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The  following  listing  summarizes  some  of  the  model  outputs 
available  either  directly  (D),  through  simple  calculations  (C),  or 
empirical  models  (M). 

1.  Meteorology:   stratus  clouds  (D),  fog  (D), 
wind-temperature-devpoint-RH  (D),  mixed  layer  height  (D). 

2.  Optical  propagation:   Cn2(M),  aerosol  extinction  (M), 
cloud  extinction  (M),  background  radiance  (C,M) 

3.  EM  Propagation:   M  profiles  (C),  inversion  duct  strength  (C), 
evaporation  duct  height  (M) 

4.  Atmospheric  dispersion:   mixed  layer  height  (D),  stability 
class  (C),  wind  variances  (M) 

5.  Remote  sensing:   temperature-humidity  profiles  (D),  aerosol 
profiles  (M),  clouds  (D) 

6.  Oceanographic:  surface  fluxes  (D),  radiation  budget  - 
shortwave  and  longwave  (D,M) 
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Figure  7.   Schematic  representation  of  the  atmospheric  aerosol 
distribtuion  in  the- marine  regime.   The  sea  salt 
(locally  generated  by  whitecaps)  aerosols  (Vs)  are 
confined  to  the  mixed  layer  while  the  global  backgro 
continental  aerosols  (Vc)  are  present  above  and  belo 
the  inversion. 
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VII.   VERIFICATION 

Since  we  have  not  outlined  a  specific  integrated  type 
boundary  layer  model,  but  have  discussed  in  more  general  terms  the 
various  components  available,  the  question  of  model  verification 
becomes  somewhat  nebulous.   It  is  possible  to  list  examples  from 
the  literature  where  specific  models  have  been  tested  against 
actual  sets  of  data.   Since  people  rarely  publish  negative 
results,  these  papers  generally  can  be  described  as  "verification" 
subject  to  certain  caveats  about  the  nature  of  the  agreement.   For 
example,  if  the  divergence  is  not  available,  one  can  find  the 
divergence  value  that  gives  the  best  agreement  with  the  changes  of 
h  over  the  period  of  interest  and  then  compare  the  model 
predictions  of  the  well-mixed  6  and  q  value  to  the  measurements. 
This  is  more  a  test  of  the  self  consistency  of  the  model  than  a 
test  of  the  model's  predictive  capability.   If  the  synoptic  scale 
forcing  terms  are  accurate  however,  the  verification  is  an 
examination  of  the  integrated  mixed  layer  assumption,  and  the 
physical  parameter izations  used.   In  view  of  this,  we  have 
arbitrarily  designated  a  set  of  verification  categories  that  are 
listed  in  order  from  the  ideal  (1)  to  the  less  than  ideal. 

1)  The  model  was  intitialized  with  data  available  at  the 
start  time  and  used  to  make  a  prediction  which  was 
subsequently  verified  to  some  stated  accuracy  in  the 
atmosphere . 

2)  A  complete  data  set  was  obtained  that  allowed  an  initial- 
ization and  specification  of  synoptic  scale  forcing.   The 
model  was  found  to  accurately  reproduce  the  atmospheric 
measurements . 
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3)  The  model  was  verified  on  a  set  of  data  obtained  in 
laboratory  experiments  where  divergence  and  advection 
(presumably)  are  negligible  or  known. 

4)  The  model  was  verified  with  an  adjustable  synoptic 
scale  forcing  parameter. 

5)  The  model  was  verified  with  an  adjustable  physical 
parameterization. 

A  summary  of  some  recently  published  papers  is  given  in  Table 
I.   The  results  might  optimistically  be  termed  "encouraging".   Not 
suprisingly,  there  is  no  Type  1  verification.   Only  one  of  the 
Type  2  verifications  was  made  for  the  marine  atmospheric  regime 
and  that  particular  example  did  not  examine  the  time  evolutions 
but  merely  verified  that  the  model  gave  an  approximate  "steady- 
state"  in  agreement  with  the  observations.   Obviously  a  definitive 
and  comprehensive  verification  of  these  models  has  not  been  done. 
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Figure  1.   Block  diagram  of  coupled  large  scale  and  boundary  layer 
models  for  the  atmosphere  and  ocean. 

Figure  2.   Schematic  representation  of  the  two  layer  model 
structure  for  clear  sky  conserved  atmospheric 
variables-water  vapor  mixing  ratio,  Q/P,  and  virtual 
potential  temperature,  9y.   The  height  of  the 
well-mixed  layer  is  h   capped  by  an  inversion  of 
thickness  cth. 

Figure  3.   Cloud  droplet  spectral  parameters  from  mid-Atlantic 

stratus  (Ref.13).   The  dispersion  is  the  ratio  of  the 
standard  deviation  to  the  mean  radius  of  the  size 
spectrum. 

Figure  4.   Net  shortwave  (darkened  boxes)  and  longwave  (clean 

boxes)  heating  rates  (flux  divergences)  measured  from 
an  aircraft  are  shown  on  the  left.   The  right  side 
depicts  the  total  radiation  heating  rate  profile.   The 
insert  compares  these  measurements  with  two  recent 
models  (Ref.13). 

Figure  5.   The  entrainment  rate,  We ,  as  a  function  of  total 

integrated  cloud  liquid  water  content,  W  (Eq.  16a). 
The  influence  of  solar  radiation  is  shown  for  three 
different  solar  zenith  angles  (Ref.  7).   Note  8  =  90 ° 
correspond  to  zero  solar  flux. 

Figure  6.   Thermal  advection  heating  rate,  -U  •  79,  as  a 

function  of  altitude.   The  data  are  obtained  from  four 
rawindsonde  stations  in  the  southern  California  coastal 
region.   The  solid  line  is  calculated  using  the  thermal 
wind  equation  (Eq.32)  and  the  dashed  line  is  from  the 
gradient  method  (Eq.26).   The  close  agreement  of  the 
two  methods  is  not  typical. 

Figure  7.   Schematic  representation  of  the  atmospheric  aerosol 
distribtuion  in  the  marine  regime.   The  sea  salt 
(locally  generated  by  whitecaps)  aerosols  (Vs)  are 
confined  to  the  mixed  layer  while  the  global  background 
continental  aerosols  (Vc)  are  present  above  and  below 
the  inversion. 
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Atmospneric  Marine  Boundary  Layer  Predictions  for  Weapons/Systems 


Abstract 

Weapons/properties  of  the  lower  atmosphere  are  very  important  in  the  opera- 
tion of  many  systems.  Of  particular  importance  is  the  affect  on  electromagnetic 
propagation  for  the  full  range  from  optical  to  radio  wavelengths.  Small  scale 
index  of  refraction  fluctuations  and  aerosols  severely  affect  optical  propaga- 
tion. Large  scale  changes  in  the  index  of  refraction  cause  bending  of  RF  and 
microwaves,  leading  to  ducting,  fading,  and  a  host  of  other  phenomena.  Experi- 
mental assessment  of  atmospheric  properties  can  only  be  done  occasionally  and  on 
a  widely  separated  basis.  Thus,  it  is  critical  for  operational  planning  to  have 
available  a  model  which  can  extend  routine  atmospheric  measurements  (radiosonde) 
in  both  time  and  space.  At  NPS  extensive  work  is  underway  on  modeling  the  marine 
atmospheric  boundary  layer.  We  sketch  here  the  basics  of  the  model  and  describe 
how  it  can  be  applied  to  a  wide  range  of  military  problems.  We  illustrate  the 
features  being  described  and  the  present  prediction  capabilities. 
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1.0  Introduction 

Modern  warfare  has  become  critically  dependent  on  the  entire  elec- 
tromagnetic spectrum  for  command  and  control  communications,  for  weapon 
guidance,  for  electronic  warfare  support  and  for  countermeasures.  Tacti- 
cally essential  systems  are  highly  affected  by  the  environment  even  when 
conditions  are  not  severe  in  the  historical  sense.  Enhancement  or  degra- 
dation of  EM/EO  system  operational  performances  has  become  a  primary  con- 
cern of  task  force  commanders.  As  such,  the  deployment  of  resources  and 
the  modification  of  tactics  based  on  environmental  factors  in  EM/EO  pro- 
pagation will,  to  a  very  large  extent,  determine  the  effectiveness  of 
sensor,  weapon  and  communication  systems. 

Environmental  effects  on  EM/EO  system  performance  can  be  grouped  in 
the  general  categories  of  refraction  (EM)  ,  wave  front:  distortion  (EO)  , 
extinction  (EO) ,  and  particulate  dispersion  (EM/EO).  Environmental  fac- 
tors contributing  to  these  effects  are  the  vertical  gradients  of  tempera- 
ture and  humidity  for  refraction,  small  scale  inhomogeneities  (turbu- 
lence) of  the  index  of  refraction  for  optical  wave  front  distortion,  con- 
centrations of  water  vapor  and  aerosols  for  extinction  (absorption  and 
scattering)  and  turbulent  transport  for  dispersion.  All  of  these  factors 
are  multi-variable  in  terms  of  the  dependence  on  routinely  measured  and 
predicted  meteorological  variables;  pressure,  wind,  temperature  and  mois- 
ture. The  only  factors  for  which  accuracies  of  existing  and  foreseen 
measurements  approach  direct  descriptions  of  the  effect  is  that  of  gra- 
dients of  index  of  refraction  and  absorption.  Turbulence  and  aerosol  de- 
scriptions will  have  to  be  obtained  by  indirect  methods  due  to  measure- 
ment complexities  which  preclude  direct  measurement  on  operational  ships. 

We  are  interested  in  describing  all  these  factors  for  the  atmo- 
spheric region  extending  from  the  surface  to  .5  -  2  km  above  the  surface, 
the  marine  atmospheric  boundary  layer  (MABL) .  The  MABL  is  cooler  and 
more  moist  than  the  overlying  air  and  is  capped  by  a  layer  (inversion) , 
50-100  meters  thick,  in  which  temperature  increases  and  humidity  de- 
creases with  height.  Critical  values  of  the  index  of  refraction  gradi- 
ents leading  to  EM  propagation  anomalies  may  exist  within  the  shallow 
capping  inversion  because  of  the  humidity  and  temperature  gradients.  The 
entire  affected  region  (duct)  extends  below  the  inversion  and  determining 
its  lower  boundary  is  essential.  High  turbulence  occurs  in  the  inversion 
as  well  as  in  the  surface  layer  and  can  cause  degradation  of  horizontal 
and  slant  path  optical  propagation.  Inversion  turbulence  intensities  are 
usually  1  to  2  orders  of  magnitude  greater  than  those  below  and  2  to  3 
orders  of  magnitude  greater  than  those  above.  Tactically  significant  ex- 
tinction of  infrared  energy  is  due  to  water  vapor  absorption  and  marine 
aerosol  scattering  in  the  presence  of  high  humidity  and,  thus,  is  usually 
restricted  to  the  MABL. 
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All  of  the  above  features  depend  critically  on  the  detailed  vertical 
structure  of  the  atmosphere;  so  regional  climatologies  are  not  useful  for 
operational  predictions.  This  is  because  the  significant  vertical  struc- 
tures are  lost  in  averaging,  and  also  because  most  historical  measure- 
ments do  not  relate  to  specific  EM/EO  requirements.  Furthermore,  the 
strengths  of  the  gradients  and  their  levels,  and  hence,  the  affected  re- 
gions cannot  be  explained  on  the  basis  of  large  scale  atmospheric  flows 
since  they  are  controlled  by  near  surface  dynamic  processes  as  well  as 
larger  scale  features. 

We  believe  a  'gap'  has  existed  in  past  efforts  to  characterize  tact- 
ical environmental  conditions.  The  gap  was  between  two  extreme  approaches 
of  relating  conditions  to  1)  near  surface  observations,  and  2)  to  larger 
scale  predicted  synoptic  patterns.  Clearly,  to  assess  the  above  features, 
local  measurement  is  desirable.  Measurements  (radiosondes)  are  made  in- 
frequently and,  as  time  progresses,  the  initial  point  measurement  becomes 
less  applicable  and  a  predictive  scheme  is  needed.  One  must  consider  a 
transition  to  climatology,  large  scale  numerical  analysis  predictions,  or 
dynamic  models  based  on  the  initial  soundings  which  are  available  at  the 
operational  location.  Climatology  has  already  been  addressed  and  discoun- 
ted. Ruggles  (1975)  has  examined  the  capabilities  of  large  scale  numeri- 
cal procedures  at  remote  facilities  and  argued  convincingly  that  they  are 
not  sufficient. 

Needed  are  characterizations  of  coupled  local  oceanic  and  atmo- 
spheric mixed  layer  features  for  time  scales  from  12  to  18  hours  and  for 
spatial  scales  from  50  to  300  km.  Significant  changes  occur  in  both  mixed 
layers  over  these  temporal  and  spatial  scales.  Predicting  these  changes 
is  now  possible  using  available  measurements  and  reasonable  physical 
models  which  have  been  recently  formulated. 

This  paper  will  describe  results  from  observational  and  model  evalu- 
ation/reformulation efforts  which  provide  evidence  that  tactically  rele- 
vant forecasts  are  possible  over  the  required  time  and  space  scales.  The 
MABL  will  be  considered  but,  as  will  be  seen,  the  MABL  model  must  eventu- 
ally be  coupled  with  an  Ocean  Boundary  Layer  (OBL)  model  to  adequately 
describe  the  changes  in  both  layers. 

The  long  term  objective  of  work  in  these  areas  is  to  make  reliable 
6-18  hour  forecasts  of  properties  within  the  MABL  and  CBL.  The  near  term 
objectives  are  1)  to  verify  and  reformulate  existing  models  for  the  re- 
sponsible physical  processes  and  2)  to  develop  measurement  procedures 
which  are  suitable  for  routine  observations  from  a  single  ship  or  air- 
craft. It  will  be  necessary  for  the  finalized  model  to  predict  such 
properties  as  the  height  of  the  inversion,  the  strength  of  the  gradients 
and  fluxes  at  the  inversion  and  at  the  surface,  and  the  values  of  well 
mixed  properties  in  the  region  between  the  surface  and  the  inversion. 
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2.0  Background 

This  section  will  provide  a  general  outline  of  the  required  descrip- 
tions and  the  models  which  can  be  used  to  predict  them.  The  descriptions 
and  models  are  delineated  in  Figure  1,  where  the  descriptions  are  en- 
closed by  ellipses  and  models  and  measurements  are  enclosed  by  rect- 
angles. The  meteorological  descriptions  and  models  are  enclosed  by  a 
dashed  outline.  Tactical  descriptions  (forecasts)  to  which  the  improved 
understandings  would  apply  appear  on  the  extreme  right  hand  side.  Single 
station  assessment  constraints  associated  with  the  local  atmospheric  and 
ocean  models  are  timely  and  compatible  with  the  concept  of  the  Tactical 
Environmental  Support  System  (TESS)  as  recently  described. by  NEPRF  scien- 
tific personnel  (NEPRF  METRO  Report,  July  1981). 

As  stated,  we  will  only  consider  the  local  atmospheric  mixed  layer 
model  and  associated  descriptions.  The  local  oceanic  descriptions  would 
be  made  by  an  oceanic  mixed-layer  model  (e.g.,  Garwood,  1977)  and  the  re- 
gional descriptions  would  be  made  by  a  three  dimensional  atmospheric 
(e.g.  NORAPS)  and  a  regional  oceanic  (e.g.  TOPS)  numerical  prediction 
model.  The  ocean  mixed  layer  model  and  the  regional  numerical  prediction 
schemes  are  objectives  of  other  Navy  sponsored  basic  and  exploratory 
research  efforts. 

2. 1  Meteorological  descriptions  and  models 

From  a  local  assessment  perspective,  let  us  consider  an  idealization 
of  the  oceanic-atmospheric  system.  The  sea-air  interface  is  bordered  by 
oceanic  and  atmospheric  turbulent  mixed  layers  which  effectively  insulate 
the  bulk  ocean  and  atmospheric  regions.  The  primary  sources  of  the  tur- 
bulence within  the  layers  are  the  velocity  (current)  and  buoyancy  (den- 
sity) gradients  at  the  interface.  Even  under  conditions  where  the  water 
is  slightly  cooler  than  the  air,  buoyancy  forced  velocity  fluctuations 
within  the  layer  can  be  quite  large  and  mix  the  entire  MABL  from  the  sur- 
face to  the  inversion.  The  large  vertical  mixing  yields  constant  (well- 
mixed)  wind,  temperature  and  humidity  profiles  above  the  surface.  At  the 
top  of  the  atmospheric  mixed  layer  there  is  a  thin  transition  region  (in- 
version) .  Examples  of  observed  well-mixed  profiles  and  the  inversion  are 
shown  in  Figure  2(a). 

The  atmospheric  mixed  layer  interacts  with  the  free  atmosphere  at 
the  inversion  by  means  of  turbulence  forced  entrainment.   Entrainment 
brings  dry,  warm  air  into  the  mixed  layer  and  also  increases  the  surface 
winds  if  there  are  higher  winds  aloft.   Stratus  clouds  would  form  within 
the  layer  if  entrainment  caused  the  moist  mixed  layer  to  extend  above  the 
condensation  level. 

Physical  processes  within  the  mixed  layer  are  also  controlled  or 
forced  (over  a  time  scale ,of  -24-36  hours)  by  the  large  scale,  non-sta- 
tionary, synoptic  flows.  Synoptic  scale  mean  vertical  motion  and  advec- 
tion  will  have  to  be  prescribed  or  predicted  ever  the  forecast  period. 
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The  previously  inferred  well-mixed  nature  of  the  convective  MABL  has 
implications  for  the  vertical  distributions  of  mean  values  and  the  verti- 
cal fluxes  (transports)  of  wind,  temperature  and  humidity.  One  implica- 
tion already  illustrated  is  that  properties  which  are  conserved  during 
mixing  can  be  treated  as  being  constant  with  height  within  the  MABL. 
These  parameters  are  the  specific  humidity  and  potential  temperature  for 
a  clear  MABL  and  equivalent  potential  temperature  and  total  (vapor  plus 
liquid)  specific  humidity  for  a  cloudy  MABL.  A  second  implication  is 
that  vertical  fluxes  of  the  well-mixed  parameters  decrease  linearly  with 
height. 

These  implications  enable  predictions  of  MABL  evolution  to  be  based 
solely  on  fluxes  at  the  upper  (inversion)  and  lower  (surface)  boundaries 
and  the  large  scale  subsidence  and  advection.  They  form  the  basis  of 
recent  model  formulations  by  Deardorff  (1976)  and  Stage  and  Businger 
(1981).  Fluxes  of  both  boundaries  are  due  to  buoyant  and  mechanical 
generated  turbulence  for  the  clear  case  with  the  addition  of  cloud  top 
radiative  fluxes  for  the  cloudy  case.  The  linear  height  variations  of 
the  fluxes  allow  one  to  relate  buoyant  fluxes  at  the  inversion  to  the 
more  readily  estimated  surface  fluxes.  Also,  cloud  top  radiative  fluxes 
can  be  estimated  on  the  basis  of  general  cloud  features.  Approaches  exist 
for  estimating  synoptic  scale  forcing  from  single  station  measurements 
but  further  efforts  are  required  to  achieve  the  accuracies  required  in 
MABL  predictions. 

In  general,  existing  models  are  quite  good  for  the  clear  sky  MABL 
and  fair  to  good  for  the  cloudy  MABL.  Considerable  effort  is  now  being 
directed  to  improving  the  models  for  the  cloudy  MABL.  The  existence  of 
cloud  layers  in  the  MABL  is  very  sensitive  to  both  the  sea  surface  tem- 
perature and  the  entrainment  rate  and,  in  turn,  clouds  have  a  profound 
effect  on  the  short  and  long  wave  radiation  budget  at  the  surface. 
Fairall  et  al  (1982)  provide  our  position  discussion  on  the  approach  and 
the  status  of  abilities  to  predict  the  evolution  of  the  MABL. 

2. 2  Tactical  descriptions 

Figures  2b-d  illustrate  tactical  descriptions  discussed  in  the  In- 
troduction and  appearing  in  the  right  hand  column  of  Figure  1.  Meteoro- 
logical processes  and  features  relevant  to  these  tactical  descriptions 
are  those  that  were  discussed  above  (Figure  2a  and  column  4  of  Figure  1) . 
Several  meteorological  factors  affect  each  tactical  description  as  indi- 
cated by  the  numbers,  1-4,  in  Figure  1.  A  few  of  the  significant  proper- 
ties and  values  in  the  tactical  descriptions  are  described  in  the 
following  paragraphs. 

The  M  profile  in  Figure  2b  describes  the  refracted  EM  ray  radius  of 
curvature,  relative  to  the  earth's  radius  of  curvature.   If  M  decreases 
with  height,  a  trapping  layer  is  formed  because  the  EM  ray  will  bend 
downward  relative  to  the  earth.  This  causes  the  formation  of  a  duct.  The 
upper  boundary  of  the  duct  is  at  the  M  minimum;  the  lower  boundary  is  at 
the  height  where  this  same  value  of  M  occurs  below  the  trapping  layer,  as 
shown  in  Figure  2b.  The  M  profile  is  determined  by  pressure,  temperature 
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and  humidity  profiles,  so  values  of  well  mixed  temperature  and  humidity 
determine  the  lower  boundary  of  the  duct.  An  increase  of  mixed  layer  M 
value  with  time  could  cause  an  elevated  to  become  a  surface  based  duct  as 
shown  in  Figure  2b.  This  could  occur  if  the  mixed  layer  became  warmer 
and  drier  due  to  entrainment  of  overlying  air.   (Warming  and  drying  of 
the  mixed  layer  by  entrainment  would  also  increase  the  height  of  the 
evaporation  duct,  a  ducting  layer  immediately  adjacent  to  the  surface.) 
Examples  of  extended  ranges  and  holes  occurring  with  a  duct  appear  in 
Figure  3  (Hitney,  1979). 

2 

The  C.  profile  in  Figure  2c  shows  the  vertical  variation  of 

turbulence  which  would  affect  optical  wave  front  distortion.  -We  see  that 

2 

C,  is  largest  near  the  surface  and  in  the  inversion.  Values  in  these 

regions  are  1-2  orders  of  magnitude  larger  than  those  in  the  mixed  layer 

2  -14  -2/3 

and  above  the  inversion.  C,  values  near  10    m    are  representative 

for  these  two  regions.  The  importance  of  this  value  is  illustrated  in 

Figure  4  with  a  simulated  image  of  a  remotely  piloted  vehicle  as  viewed 

2 

through  C.     region 

and  Walter,  1978). 


through  C.  regions  of  0  (no  turbulence)  ,  10   ,  and  10    m~    (Kearns 


The  extinction  coefficient,  8,  profile  is  shown  in  Figure  2d.  3  is 
inversely  proportional  to  the  range  of  IR  systems  and  depends  on  the 
absolute  humidity  and  the  aerosol  size  distribution.  Aerosol  size  dis- 
tributions depend  on  the  generation  rate  of  sea  salt  particles  and  on  the 
relative  humidity  through  the  growth  factor.  Therefore,  IR  extinction  is 
very  dependent  on  the  entrainment  of  relatively  clean,  dry  and  warm  air 
from  above  the  inversion  and  on  the  surface  wind  generation  of  the  sea 
salt  particles.  As  illustrated,  the  extinction  coefficient  is  large 
within  the  mixed  layer  and  increases  with  height  due  to  relative  humidity 
increases.  The  increase  with  height  is  important  in  slant  path  range 
considerations. 

3. 0  MABL  observations  and  predictions 

A  recently  obtained  data  set  will  be  used  to  demonstrate  the  nature 
of  changes  in  the  MABL  and  their  prediction.  The  demonstration  will 
show: 

a)  the  status  of  available  data  sets, 

b)  that  the  assessment  of  relative  roles  of  dynamic  processes  in 
the  MABL  can  be  based  on  a  relatively  simple  physical  model 
which  includes  specifiable  inputs,  and 
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c)  that,  at  the  present  time,  we  have  only  initial  evidence  that 
changes  can  be  predicted.  The  rigorous  specification  or 
description  will  require  further  interpretive  efforts  and, 
perhaps,  improved  formulation  of  the  separate  models. 

3.1  MABL  observations 

The  data  were  obtained  during  the  Cooperative  Experiment  on  West 
Coast  Oceanography  and  Meteorology  (CEWCCM-78)  conducted  west  of  San 
Nicolas  Island,  CA  during  May  of  1978.  Observations  of  the  oceanic  and 
atmospheric  mixed  layers  were  made  from  the  R/V  ACANIA;  radiosonde  obser- 
vations were  also  made  at  surrounding  shore  stations.  The  data  to  be 
shown  are  from  a  48-hour  period,  5/19/1200  to  5/21/1200  PST.  The  R/V 
ACANIA  was  cruising  slowly  (2-3  knots)  into  the  wind  and  returned  to  an 
initial  point  approximately  every  12  hours.  The  general  location  of  the 
R/V  ACANIA  during  the  5/19  to  5/22  periods  and  locations  of  surrounding 
shoreline  radiosonde  sites  appear  in  Figure  5. 

The  period  was  one  of  steady  onshore  flow  caused  by  the  combined  ef- 
fects of  intensification  of  the  Eastern  Pacific  High  and  the  persistence 
of  the  Mexican  thermal  low.  The  only  apparent  change  in  synoptic  scale 
forcing  was  an  increase  in  the  offshore  pressure  gradient.  Advection  in 
the  atmosphere  was  moderate.  However,  MABL  evolutions  were  primarily 
determined  by  subsidence,  surface  fluxes,  and  entrainment  at  the 
inversion.  Satellite  imagery  showed  increasing  uniform  stratus  coverage 
(thin  to  heavy)  during  the  period  with  a  cellular  (broken)  coverage 
occurring  late  on  the  21st. 

Acoustic  sounder  returns,  mean  surface  layer  parameters,  and  sea 
surface  temperature  during  the  19-21  May  period  are  shown  in.  Figures 
6b-d,  all  measured  from  the  R/V  ACANIA.  Potential  temperature  composite 
profiles  from  shipboard  and  shore  station  radiosonde  and  temperature  pro- 
files from  ship  deployed  X3T's  appear  in  Figures  6a  and  6e,  respectively. 
Although  these  descriptions  of  MABL  and  CBL  changes  were  obtained  from  an 
instrumented  research  vessel,  all  can  be  obtained  from  operational  ships. 
This  includes  the  acoustic  sounder  record. 

Changes  of  atmospheric  features  which  would  have  been  tactically 
significant  in  view  of  the  time  scale  and  their  magnitudes  were: 

a)  The  MABL  depth  increased  from  250  m  to  750  m  over  the  48  hour 
period.  The  changes  occurred  in  relatively  short  intervals, 
from  20/00  to  20/04,  and  from  21/00  to  21/08.  The  ABL  depth 
remained  near  constant  during  the  intervening  12  to  18  hour 
periods. 

b)  The  surface  layer  temperature  increase  from  20/00  to  20/15  is 
indicative  of  entrainment  of  overlying  warm  air.  Entrainment 
would  also  dilute  the  mixed  layer  and  decrease  the  relative 
humidity.  The  entrainment  was  presumably  also  a  factor  in  the 
wind  speed  increase  over  the  same  period,  which  caused  the  warm 

A-8 


UNCLASSIFIED 

shallow  ocean  layer  to  be  destroyed,  lowering  the  sea  surface 
temperature  after  20/1800  (Figure  6d) . 

In  conjunction  with  changes  in  f-PiBL  structures  and  parameter  values 
depicted  in  Figure  6,  changes  occurred  which  are  important  to  operational 
systems.  Some  of  these  changes  can  be  determined  directly  from  parameters 
measured  in  this  scientific  observational  effort.  M  profile  and  cloud 
base  heights  (condensation  levels)  can  be  calculated  from  the  measured 
humidity  and  temperature  profiles.  Surface  layer  extinction  coefficients 
have  to  be  determined  from  the  measured  humidity  and  aerosol  size  distri- 
butions. Our  approach  for  determining  extinction  coefficients  from  aero- 
sol data  has  been  described  by  Schacher  et  al  (1981).  Surface  .layer  C, 
values  and  evaporation  duct  heights  can  be  determined  with  considerable 
certainty  from  measured  surface  layer  wind,  temperature,  humidity  and 
surface  temperature  values  using  expressions  described  by  Davidson  et  al 
(1981)  and  Fairall  et  al  (1978),  respectively.  Extinction  and  C.  pro- 
files have  to  be  based  on  more  recent  and,  hence,  less  substantiated  ex- 
pressions such  as  those  by  Wells  et  al  (1977)  and  Wyngaard  and  Lemone 
(1980)  . 

Results  from  these  determinations  appear  in  Figures  7  and  8  along 
with  the  composite  potential  temperature,  3,  and  water  vapor  mixing 
ratio,  q,  profiles  (Figure  7a)  and  the  inversion  height  (Figure  8a,  rect- 
angles). The  observed  results  are  the  solid  lines  for  profiles  (Figure 
7)  and  the  X's  for  observed  (Figure  8).  The  observed  inversion  heights 
(Figure  8a)  (the  rectangles)  are  shown  for  all  stations  used  in  the 
composites.  The  observed  results  indicate  that  for  the  observation 
period: 

a)  The  EM  duct  associated  with  the  inversion  gradients. evolved  from 
being  surface  based  to  being  elevated  (Figure  7b)  and  the  eva- 
poration duct  height  ranged  from  below  5  m  to  above  3  m  (Figure 
8b), 

2  -15 

b)  Optical  turbulence,  CL,  ,   in  the  surface  layer  varied  from  10 

to  10         m  over  the  period    (Figure  8c)    and  values  in  the 

inversion  increased  from  10~      to  10        m  (Figure  8c) , 

c)  Surface  layer  aerosol  extinction  in  the  8-12  urn  IR  region  in- 
creased during  the  first  6  hours    (5/19  1200-1800)    because  the 
relative  humidity  increased.      It  decreased  during  the  last  18 
hours   (5/20/1800  to  5/21/1200)    because  of  relative  humidity  and 
wind  decreases   (Figure  8d) , 

d)  stratus  cloud  persisted  through  the  period.     The  base  lifted 
from  about  100  m  at  5/19/1300  to  600* meters  at  5/21/1200.     This 
is  evident  in  the  extinction  profile    (Figure  7d)  . 
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3.2  Model  predictions 

Daring  the  past  five  years,  we  performed  investigations  on  surface 
layer  and  mixed- layer  scaling  of  small  scale  turbulence  and  aerosol  pro- 
perties (Davidson  et  al,  1978;  Fairall  et  al,  1979;  Schacher  et  al, 
1981) .  The  investigations  led  to  the  application  of  a  time  dependent 
MABL  model  based  on  entrainment  energetics  and  on  cloud  radiative  trans- 
fers suggested  by  Deardorff  (1976) .  The  model  requires  as  input  an  ini- 
tial atmospheric  sounding  (radiosonde) ,  the  mean  winds  at  a  level  within 
the  surface  layer  (10-30  meters)  and  the  surface  temperature.  Well-mixed 
temperature  and  humidity  are  predicted  so  the  surface  wind  and  wind  shear 
at  an  inversion  are  the  only  local  atmospheric  variables  which  have  to  be 
prescribed.  The  larger  scale  subsidence  and  advection  must  also  be 
prescribed  for  the  forecast  period. 

The  steps  in  the  prediction  computation  are  shown  in  Figure  9,  where 
it  is  noted  that  procedures  are  the  same  for  clear  and  cloudy  cases  ex- 
cept for  entrainment  computation  and  estimating  cloud  top  cooling. 
Fairall  et  al  (1982)  include  comprehensive  discussions  of  how  responsible 
physical  processes  are  treated  in  the  model  and  how  available  data  can  be 
used  to  estimate  the  processes.  Because  of  the  simplified  physical  model 
the  computations  do  not  require  numerical  integrations  on  a  vertical  grid 
so  the  computer  storage  requirements  can  be  satisfied  by  available  ship- 
board microcomputers.  In  fact,  we  are  examining  the  use  of  a  hand-held 
calculator  for  the  procedures. 

Predictions  of  MABL  changes  were  made  with  the  simplified  model  for 
the  48  hour  period  (5/19/1200  to  5/21/1200)  corresponding  to  observed  re- 
sults in  Figures  6a-e.  The  predictions  were  made  for  two  separate  24  hour 
periods  starting  at  19/1200  and  at  20/1200.   Initial  profiles  for  the 
19/1200  start  time  was  that  obtained  at  San  Clemente  Island  (SCI  in 
Figure  5)  at  approximately  1000.   Initial  profiles  for  the  20/1200  start 
time  were  averages  of  the  composites  at  20/0500  and  20/1700.  The  sea  sur- 
face temperatures  were  those  at  each  start  time  as  shown  in  Figure  6d  and 
were  15  and  14.5  C.  The  surface  wind  observed, during  the  forecast  period 
was  specified  to  increase  linearly  from  5ms  .at  19/1200  to  10  m  s   at 
20/2100  and  then  to  decrease  linearly  to  7  m  s   at  21/1200. 

Larger  scale  synoptic  forcing  (subsidence  and  thermal  advection) , 
stated  previously  as  requiring  more  study  if  it  is  to  be  satisfactorily 
estimated  from  single  station  assessments,  was  estimated  from  temperature 
changes  above  the  inversion  (subsidence)  and  from  the  thermal  wind 
(thermal  advection) .   It  was  set  to  zero  during  the  first  period.  Subsi- 
dence for  the  second  period,  starting  at  20/1200,  was  adjusted  on  the 
basis  of  that  value  required  for  agreement  between  predicted  and  observed 
mixed  layer  depths  during  the  first  period.  Hence,  the  model  was  also 
used  to  estimate  the  "most  probable"  synoptic  scale  subsidence. 

The  predictions  are  illustrated  in  Figures  7  and  8  along  with  the 
"observed"  results.  The  predictions  are  the  dashed  lines  for  both  the 
profiles  (Figure  7)  and  inversion  height  and  surface  layer,  10  m,  values 
(Figure  8) .  The  agreement  near  20/1200  in  Figure  8  occurs  because  the 
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model  was  reinitialized  at  that  point.  Significant  outcomes  from  compar- 
isons of  observed  and  predicted  results  are 

a)  The  mixed  layer  depth  (inversion  height,  Figure  8a)  was  predic- 
ted quite  well.  This  increase  during  the  later  part  of  both 
periods  was  caused  by  cloud  top  entrainment  because  the  predic- 
ted surface  layer  stability  became  stable  and,  hence,  there 
existed  no  surface  buoyancy  flux  forcing. 

b)  The  cloud  cover  was  accurately  predicted  to  persist  throughout 
the  first  period  and  for  the  first  15  hours  of  the  second 
period,  as  indicated  by  the  condensation  level  being  below  the 
inversion  in  Figure  8a.  However,  the  predicted  cloud  thickness 
was  less  than  that  which  actually  occurred.  As  indicated  in 
preceding  discussions,  the  cloud  cover  was  observed  to  break  up 
in  the  21st  and  the  prediction  after  21/0600  agrees  with  that. 

c)  The  well  mixed  temperature  and  humidity  and  the  jump  predictions 
(Figure  7a)  are  reasonable  except  for  the  20/0500  comparison 
time.  The  predicted  well  mixed  temperature  was  too  warm  (1-2  C) 
at  all  other  comparison  times.  This  led  to  the  surface  layer 
becoming  stable  when  in  actuality  it  remained  unstable  through- 
out both  periods,  surface  temperature  higher  than  10  m  tempera- 
ture. We  believe  the  cooling  associated  with  the  clouds  was  not 
sufficient.  It  will  be  seen  that  the  surface  layer  was  stable 
(too  warm)  during  the  night  hours. 

d)  The  predicted  M  profiiles  (Figure  7b)  were  very  accurate  in 
terms  of  the  surface  based  duct  evolving  to  an  elevated  duct  and 
the  lower  boundary.  It  is  interesting  that  even  when  the  pre- 
dicted well  mixed  temperature  and  humidity  were  much  too  high 

(20/0500)  the  M.  profile  verified.  This  occurred  because  their 
respective  affects  on  M  compensated. 

2 

e)  The  predicted  C.  profiles  (Figure  7c)  gave  good  values  for  the 

inversion  region  and  reasonable  values  for  the  mixed  layer,  ex- 

•  cept  for  20/0500.   It  is  noted  that  values  above  the  inversion 

—18 
were  set  to  10  ~  for  both  observed  and  predicted  profiles. 

2 
Factors  of  3  agreement  are  considered  reasonable  for  C.  .   It  is 

noted  that  inversion  region  values  should  be  based  on  the  con- 

vective  mixing  velocity,  WA  (Wyngaard  and  Lemone,  1980).  When  a 

stable  surface  layer  existed  (20/0500)  the  friction  velocity, 

U^,  was  used  instead  of  W+. 
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f)  The  total  extinction  at  8-12  urn  prediction  (Figure  7d)  was  quite 
good  outside  of  cloud  layers.  Again,  values  above  the  inversion 
were  set  to  the  same  value  for  both  observed  and  predicted  pro- 
files. The  sensitivity  of  aerosol  size  to  relative  humidity 
changes  when  the  relative  humidity  is  above  90%  is  what  leads 
this  comparison  to  give  seemingly  poor  results.  Hence,  the  good 
agreement  outside  of  clouds  is  encouraging. 

g)  The  evaporation  duct  height  (Z*)  prediction  (Figure  8b)  is  not 
very  good.  The  disagreement  occurs  because  the  predicted  sur- 
face was  less  unstable  (more  stable)  than  the  observed.  Z*  in- 
creases almost  linearly  with  increasing  stability.  The  very 
large  predicted  Z*  values  at  the  end  of  the  second  period  oc- 
curred because  it  was  stable  and  because  the  predicted  well 
mixed  humidity  decreased,  Figure  7a. 

2 

h)  The  predicted  surface  layer  (L,  values  (Figure  8c)  were  always 

less  than  those  observed  because  of  the  previously  mentioned  in- 
accurate surface  layer  stability.  The  minima  at  19/2300,  20/1800 
and  20/0500  correspond  to  the  surface  layer  going  from  unstable 
to  stable,  unstable  to  stable  and  stable  to  unstable, 
respectively. 

i)  The  predicted  surface  layer  aerosol  scattering  extinction  for 
8-12  pm  (Figure  8d)  was  very  good  during  the  second  period.  The 
observed  and  predicted  decrease  from  20/1800  to  21/1200  was 
associated  with  due  to  a  decrease  in  relative  humidity  and  wind 
speed.  The  wind  speed  used  for  the  prediction  was,  of  course, 
prescribed  on  the  basis  of  observed  results.  Hence,  the  role  of 
relative  humidity  on  the  extinction  is  that  being  evaluated. 

4.0  Conclusions 

The  importance  of  and  an  approach  for  predicting  evolutions  of  the 
marine  atmospheric  boundary  layer  have  been  described.  Any  change  in  the 
mixed  layer  depth  or  well-mixed  properties  can  lead  to  significant 
changes  in  the  weapons/systems  tactical  environment.  These  include  EM 
duct  regions,  optical  turbulence  and  optical  extinction,  both  within  and 
at  the  boundaries  of  the  inversion  capped  mixed  layer. 

Although  the  tactical  descriptions  are  multi-variable  and  require 
detailed  vertical  descriptions,  a  simplified  model  based  on  routinely 
observed  data  appears  to  be  quite  useful  for  predicting  changes  ever 
12-18  hours.  For  the  cloud  topped  period  considered,  the  refraction, 
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optical  turbulence  and  optical  extinction  profiles  were  predicted  quite 
well.  Cloud  coverge  changes  were  also  predicted  quite  well.  The  primary 
differences  between  observed  and  prediction  results  were  for  the  surface 
layer,  the  evaporation  duct  height  and  the  optical  turbulence.  The 
differences  were  significant  in  terms  of  tactical  effects.  The  predicted 
evaporation  duct  was  much  too  high  and  the  predicted  optical  turbulence 
was  much  too  low.  The  reason  in  both  cases  was  that  the  predicted  air 
temperature  was  too  high  (2-3  C)  which  resulted  in  a  neutral  to  stable 
condition.  We  believe  the  cause  for  this  was  the  improper  specification 
of  cloud  cooling  within  the  mixed  layer. 
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Fig.  2.  Simplified  meteorological  and  tactical  features  of  the  MAJ3L  and 
overlying  region  including  (a)  humidity  mixing  ratio,  Q/p,  and 
potential  temperature,  e,  profiles  and  sub-regions,  (b)  modified 

index  of  refraction,  M,  profiles,  (c)  optical  turbulence  inten- 

2 
sity,  Cjj  ,  profiles,  (d)  total  extinction  coefficient,  8 

(absorption  and  scattering) ,  profile. 
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Fig.   3.     Examples  of  extended  ranges  and  holes  for  EM  ducts  for   (a)    surface- 
based  duct  and  a  shipboard  air-search  radar  and   (b)   elevated  duct  and 
an  airborne  early-warning  radar   (from  Hitney,   1979). 
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Fig.    4.     Computer  simulated  image  remotely  piloted  vehicle    (HPV)    as  viewed  with 

8-12     m  IR  sensor    (range  2  km,  visibility  15  km,   threshold  4.5  nW) 

2 
through  atmosphere  with  C^     values  of    (a)   0    (no  turbulence),    (b)    3. 7  x 

10~15  m"2/3  and    (c)    1. 0  x  10~14  m~2/3    (from  Kearns  and  Walter,   1978). 
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Fig.  5.  Tract  of  R/V  Acania  and  location  of  shoreline  radiosonde  sites  during 
CEWCCM-78.  General  location  of  R/V  Acania  during  19-21  May  period  is 
indicated  by  hatched  area  N-MW  of  San  Nicolas  Island  (SNI). 
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Fig.   6.     Atmospheric  and  oceanic  mixed-layer  observations  during 
CEWCCM-78    (a)    potential   temperature  composite  profiles, 
(b)    acoustic  sounder  record,    (c)    10  meter  wind  speed,    (d)    10 
meter  and  surface  temperatures  and    (e)   XBT  traces. 
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Fig.    7.      Observed/computed   and   predicted  MABL  profiles   for  5/19/1200  to 
5/21/1200  CEWCCM-78  period    including    (a)    Q    (g  m~3)    and   9(°C), 

(b)    M,    (c)    C        (m     '   )    and    (d)    total   extinction  coefficient,    8 

-1 
(km     ).     Solid   lines  correspond   to  observed  and  dashed   to  model 

prediction. 
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Fig.    8.      Observed/ computed  and  predicted  surface  layer  values   for 

5/19/1200  to  5/21/1200  CEWCCM-78  period   including    (a)    inve-sion 

height  and  lifting  condensation  level,    (b)    evaporation  duct, 

2     -2/3 
(c)    CN     m  and    (d)    aerosol   scattering  extinction  coefficient, 

3.      X's  are  observed  and  dashed  lines  are  model  predictions. 
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Meteorological  models  for  optical  properties  in  the  marine  atmospheric  boundary  layer 

C.  W.  Fairall 
BDM  Corporation,  Naval  Postgraduate  School,  Monterey,  California  93940 
K.  L.  Davidson  and  G.  E.  Schacher 
Naval  Postgraduate  School,  Monterey,  California  93940 

Abstract 

Observational  experiments  on  turbulent  intensities  and  aerosol  distributions  in  the 
marine  boundary  layer  (MBL)  have  been  performed  over  several  years.   Observations  have  been 
made  with  ship-mounted  and  airplane-mounted  sensors.   Objectives  have  been  to  relate  optical 
properties  to  meteorological  descriptions  which  utilize  scaling  laws  for  the  MBL.   The 
approach  has  been  to  incorporate  in  the  descriptions  the  surface  fluxes  of  momentum,  heat, 
and  moisture,  the  processes  at  the  inversion  and  the  profiles  within  the  intervening  convec- 
tively  mixed  layer.   We  have  found  that  optical  turbulence  parameters  (C^2  and  1q )  can 
be  readily  estimated  using  measured  mean  values  of  wind,  temperature,  and  humidity  with 
recent  bulk  formulae  to  derive  the  surface  fluxes.   These  estimates  appear  to  be  more  reli- 
able than  values  obtained  from  direct  (but  difficult  to  perform) . turbulence  measurements. 
The  model  for  obtaining  the  estimates  was  evaluated  on  the  basis  of  optical  C^    values 
with  good  agreement.   Good  comparisons  have  been  observed  between  extinction  values  obtained 
from  transmission  measurements  and  those  obtained  from  calculations  on  measured  aerosol 
distributions.   Existing  empirical  formulations  which  related  the  latter  to  wind  speed  and 
relative  humidity  appear  to  be  inadeguate  except  for  climatological  purposes.   This  is 
because  other  influences  on  equilibrium  aerosol  distribution  are  not  included. 
Reformulation  of  these  expressions  is  being  performed  to  include  the  height  of  the  inversion 
(mixing  volume)  and  surface  fluxes  (aerosol  generation  and  transport). 

Introduction 

Three  atmospheric  processes  are  primarily  responsible  for  the  degradation  of  the  trans- 
mission of  optical  images  and  electro-optical  energy:   aerosol  extinction,  molecular  extinc- 
tion, and  turbulent  distortion  (scintillation  and  beam  wander).   Optical  energy  propagating 
through  the  atmosphere  is  scattered  and  absorbed  by  aerosols  and  molecules  and  optical  wave 
fronts  are  deflected  and  distorted  by  turbulent  fluctuations  in  the  refractive  index. 
Tested  physical  models  of  these  atmospheric  effects  and  approaches  to  estimate  them  from 
routine  meteorological  data  are  required  for  applications  of  optical  and  electro-optical 
systems.   Naval  Postgraduate  School  (NPS)  investigators  have  participated  in  field  experi- 
ments designed  to  improve  and  verify  such  models  and  approaches  for  the  overwater  regime. 
These  have  been  performed  for  U.  S.  Navy  research  programs  over  the  past  seven  years. 

The  two  atmospheric  optical  parameters  of  primary  interest  are  the  total  extinction  coef- 
ficient and  the  refractive-index  structure  function  parameter,  Cj,j2.   The  total  extinc- 
tion coefficient  (S  +  a)  parameterizes  the  loss  of  optical  energy  as  it  is  scattered  out  of 
the  beam  or  absorbed  by  molecular  and  particulate  constituents  of  the  atmosphere.   Thus,  it 
has  four  components:   molecular  scattering  and  absorption  CS  =  3a  +  0S)  and  aerosol 
scattering  and  absorption  (a  =  aa  +  a$) .   The  refactive- index  structure  function  parameter 
parameterizes  the  intensity  of  refactive-index  spatial  inhomogeneities  which  distort  and 
tilt  image  wave  fronts.   Both  of  these  parameters  can  be  related  to  meteorological 
parameters:   1)  CN2  to  the  intensity  of  small  scale  turbulent  fluctuations  and  2)  total 
extinction  to  the  concentration  of  certain  gases  and  to  the  size  distribution  of  aerosols. 

One  part  of  the  experiments  dealt  with  the  compatibility  of  optical  and  meteorological 
propagation  theories.   For  this,  direct  optical  measurements  of  extinction  and  scintillation 
across  an  overwater  path  were  compared  to  values  calculated  from  aerosol  and  turbulence  data 
obtained  at  the  midpoint  of  the  optical  path  from  an  instrumented  research  vessel.   Once 
these  relationships  are  validated,  optical  properties  can  be  deduced  from  detailed  meteoro- 
logical measurements  and  models. 

A  second  part  of  the  experiments  involved  attempts  to  apply  existing  models  for  the 
marine  atmospheric  surface  layer  and  mixed  layer  to  describe  the  optically  relevant  turbu- 
lent and  aerosol  properties.   The  surface  layer  is  considered  to  be  the  region  extending 
from  the  air-sea  interface  to  heights  of  40-50  meters.   Its  mean  and  turbulent  properties 
are  closely  related  to  surface  fluxes  of  momentum,  temperature  and  moisture,  and  hence,  the 
air-surface  (bulk)  differences  of  mean  properties.   The  fluxes  themselves  are  often  assumed 
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to  be  independent  of  height  within  the  surface  layer.   The  mixed  layer  extends  from  the  top 
of  the  surface  layer  to  the  base  of  the  capping  inversion.   Properties  of  the  mixed  layer 
are  determined  by  both  the  underlying  surface  fluxes  and  entrainment  at  the  inversion.   Fig. 
1  illustrates  the  meteorological  signature  characteristic  of  the  marine  boundary  layer 
(MBL) .   A  well-mixed  property  (in  this  case,  virtual  potential  temperature  and  water  vapor 
mixing  ratio)  has  a  gradient  in  the  surface  layer,  is  constant  in  the  mixed  layer  (up  to 
height,  h)  changes  abruptly  in  the  interfacial  layer  (thickness,  h)  and  acquires  a  fairly 
constant  gradient  (lapse  rate,  r)  in  the  free  troposphere.   This  structure  prevails  over  the 
open  ocean  a  majority  of  the  time,  the  most  common  exception  being  during  periods  of  frontal 
passage.   This  structure  is  due  to  turbulent  processes  driven  by  the  air-sea  interaction. 
These  same  processes  lend  themselves  to  scaling  law  parameterizations  that  allow  the 
development  of  physical  models  of  turbulence  and  aerosol  profiles  in  the  MBL. 

Comparison  of  optical  and  meteorological  values  for  C^2 
and  aerosol  scattering  extinction  coefficient 

CN2  for  optical  wavelengths  is  related  to  the  temperature  (C^2),  water  vapor 
(Cq2)  and  temperature-water  vapor  (C-tq)  structure  function  parameters  by-'- 

CN2  =  (79  x  10"6  P/T2)2  (CT2  +  0.113  CTQ  +  3.2  x  10"3  Cq2)  (1) 

where  P  is  the  pressure  in  mb,  and  T  the  absolute  temperature. 

The  separate  components  of  total  extinction  (8  +  a )  can  be  calculated  from  meteorological 
data.   Molecular  extinction  (3)  can  be  obtained  using  existing  models,  such  as  the  LOWTRAN 
model  developed  by  the  Air  Force  Geophysics  Laboratory^,  with  observed  temperatures  and 
water  vapor  concentrations.   Aerosol  total  extinction  (a)  can  be  calculated  from  the  aerosol 
spectral  density,  N(r) ,  as 


a  =   |  it 


2  E(n,A)  N(r)  dr  (2) 


where  r  is  the  particle  radius,  E(n,A)  the  total  scattering  efficiency  it  wavelength,  \,    and 
refractive-index,  n. 

A  series  of  field  experiments  to  establish  the  credibility  of  Eq.  1  and  Eq.  2  for 
practical  determination  of  optical  properties  culminated  in  the  MAGAT  experiment  in 
Monterey.   Since  the  results  have  been  reported  for  CN2  by  Davidson  et  al3  and  for 
aerosols  by  Schacher  et  al4,  a  brief  summary  only  will  be  presented  here. 

The  optical  measurements  were  made  over  a  13  km  overwater  path  and  are  described  by 
Crittenden  et  al^.   The  meteorological  measurements  were  made  from  the  R/V  Acania 
(operated  by  the  Department  of  Oceanography  of  NPS )  which  was  stationed  near  the  midpoint  of 
the  optical  path.   The  R/V  Acania  was  equipped  with  multi-level  measurement  systems  to 
measure  mean  and  turbulent  wind,  temperature  and  humidity  and  aerosol  size  distribution 
(.09  urn  to  14  urn  radius).   Descriptions  of  the  measurement  and  analyses  systems  and 
procedures  are  given  by  Schacher  et  al  4.   Two  measurements  requiring  further  description 
for  these  camparisons  are  those  of  the  temperature  structure  function  parameter,  CT2, 
and  of  the  aerosol  size  distributions. 

Meteorological  values  for  CN2  in  this  comparison  were  determined  primarily  from 
CT2  measurements  with  minor  (less  than  15%)  adjustments  for  the  water  vapor 
contributions,  C^q  and  Cq* .      The  Ct   measurements  were  made  with  two  resistance 
wires,  2.5  urn  platinum  sensors,  separated  a  distance  of  30  cm  in  the  crosswind  direction. 
CT^  values  were  obtained  from  RMS  values  of  the  sensed  temperature  differences. 

The  aerosol  spectra  were  measured  with  optical  particle  counters  made  by  Particle 
Measurement  Systems  (PMS)  of  Boulder,  Colorado.   The  total  system  consisted  of  two  probes, 
the  classical  scattering  (CSAS)  and  the  active  scattering  (ASAS),  controlled  by  a  DAS-32 
with  computer  interfacing.   This  system  measures  aerosols  in  90  size  channels  from  0.09  ym 
to  14.0  vim  radius.   The  M(r)  spectra  were  determined  for  half-hour  periods.   Each  spectrum 
was  fit  in  L0G(N(r)),  LOG  (r)  space  with  a  seventh  order  polynomial  for  0.09  urn  <  r  <  7  ym, 
with  a  linear  fit  for  r  >  7ym.   The  extinction  was  then  calculated  using  these  fits  for 
0.03  um  <  r  <  30  ym. 

Comparisons  of  optically  measured  Cjj  with  meteorologically  measured  Ckj2 
(calculated  from  C^2   measurement  with  correction  for  probable  water  vapor  contributions) 
are  shown  in  Figure  2.   The  salt  loading  effect  increases  measured  Ct2  values  for 
general  conditions  during  these  measurements.   The  extinction  Comparisons  are  shown  in 
Fig.  3. 
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Turbulence  modeling 


At  any  given  point  in  the  atmosphere,  the  turbulence  properties  are  determined  by  the 
physical  processes  dominating  that  region  of  space.   Identifying  the  important  physical 
parameters  and  neglecting  the  unimportant  forms  the  basis  of  the  scaling  law  approach.   As 
an  extreme  example,  clearly  the  sea  surface  temperature  has  little  influence  on  the 
turbulence  at  an  altitude  of  30  km.   Therefore,  the  MBL  is  loosely  partitioned  into  regions 
in  which  certain  scaling  regimes  apply;  these  regions  are  summarized  in  Table  1.   (symbols 
will  be  defined  in  the  text). 


Table  1. 


C^  scaling  regions  of  the  MBL 


region 


height 


relevant  parameters scaling  parameters 


model 


surface 

free  convection 

transition 

interfacial 

free  troposhere 


Z<0. lh 
0.1h<Z<0.5h 
0. 5h<Z<h 
h<z<  (l+a)h 
z>(l+a)h 


surface  fluxes 
surface  fluxes,  h 
-> 

jumps,  fluxes,  h 

-> 


u*,T*,Q*,Z/L 

w*»e*,M*,z/h 

w*,we,a 


bulk 

mixed  layer 

entrainment 


Surface  layer 

Monin-Obuhkov  scaling  (MPS).   Near  the  surface,  the  height  above  the  surface,  Z,  can  be 
normalized  by  the  Monin-Obuhkov  stability  length,  L.   We  can  then  represent  the  relevant 
micrometeorological  properties  in  terms  of  height  (Z),  scaling  parameters  (T*,  Q*  and  u* ) 
and  dimensionless  functions  of  5  =  Z/L, 


CT2  =  t*2  Z-2/3  f(C) 
CQ2  =  q*2  z-2/3 


-TQ 


=  r 


TQ 


*Q*Z-2/3  Al/2  f(5) 


where  T*  and  Q*  are  the  temperature  and  humidity  scaling  parameters,  f(£)  is  a 
dimensionless  function?,  r^Q  is  the  temperature-humidity  correction  parameter  (about 
0.8)  and  A  is  a  constant  (about  0.6).   The  scaling  length,  L,  is  given  by 

L  =  (T/ieg)  u*2/ (T*  +  0.61  TQ*/P) 


(3a) 
(3b) 
(3c) 


(4) 


where  u*  is  the  friction  velocity,  <    is  Von  Karman's  constant  (0.35),  g  is  the 
acceleration  of  gravity  and  p  is  the  density  of  air.   Note  that  the  scaling  parameters  are 
related  to  the  surface  fluxes  of  momentum  (t  =  pu*2),  temperature  (Q0  =  -u*T*)  and 
water  vapor  (M0  =  -u*Q*). 

Bulk  model.   Using  Eq .  1  and  Eq .  3  one  can  calculate  C^2  a^  some  height  Z  in  the 
surface  layer,  provided  values  of  T*,  Q*  and  L  are  known.   The  most  straightforward 
method  of  obtaining  these  variables  is  to  measure  the  surface  fluxes  of  momentum,  heat  and 
water  vapor.   The  difficulty  of  direct  flux  measurements  led  -o  the  development  of  a  method 
that  utilizes  bulk  meteoroloaical  quantities  (wind  speed,  u,  temperature,  T,  and  water  vapor 
density,  0).   In  this  case,  the  scaling  parameter  for  X(X  =  u,T,Q)  is  obtained  from  the 
difference  in  X  from  the  sea  surface  (Xs)  to  some  reference  height  (usually  10  m)  in  the 
atmosphere. 


X*  =  cxl/2(X1( 


X. 


(5) 


where  cx  is  the  drag  coeficient  for  X  (typically,  cx  =  1.3  x  10  3  over  the  ocean) 
Further  details  on  use  of  the  bulk  method  can  be  found  in  Davidson  et  al  °. 


The  verification  of  the  bulk  method  for  C*t<2  ancj  Eq .  3a  is  shown  ^  in  Fig.  4a  while 
the  corresponding  verification  for  Cq2  and  Eq.  3b  is  shown^  in  Fig.  4b.   Following 
this  work,  the  MAGAT  experiment  was  performed  to  verify  the  bulk  predictions  against  direct 
optical  (scintillation)  measurements  3  (Fig.  5). 

In  this  comparison  all  but  three  of  26  pairs  of  values  are  within  a  factor  of  two  of 

perfect  agreement.   The  mean  percent  error,  taking  the  optical  values  to  be  correct,  is  33% 

which  is  extremely  good  agreement.   This  agreement  is  better  than  that  obtained  in  the 

comparison  with  Cjj2  derived  from  C-^2  measurements.  Such  a  result  is  partially 

expected  because  of  inherent  measurement  errors  in  the  latter. 


Results  of  these  comparisons  and  results  of  comparison  of  bulk  and  turoulence  derived 

parameters  in  the  above  referenced  papers  clearly  demonstrates  that  surface  layer  optical 

turbulence,  C^2,  is  very  satisfactorily  described  by  existing  bulk  formulae  and  scaling 

expressions.   In  fact,  over  the  ocean  Cvj2  can  oe  more  accurately  determined  from  bulk 

model  calucations  than  from  the  more  difficult  direct  measurements  of  small  scale 
temperature  and  humidity  fluctuations. 
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Mixed  Layer. 

Free  Convection  Layer.   In  the  lower  half  of  the  mixed  layer  free  convection  scaling 
applies  lv-*'  ^r  the  scaling  parameters  are  (where  Qov  =  QQ  +  0.61TMo/p) 

9*  -  Q0/W*  (6a) 

M*  =  M0/W*  (6b) 

W*  =  (gQovh/T)1/3  (6c) 

Note  that  some  authors  use  the  height  of  the  inversion,  Z^  =  (1  +  ct/2)h,  rather  than 
the  height  of  the  well-mixed  layer,  h.   Normally,  Zj  is  about  10%  greater  than  h. 

In  this  regime,  the  structure  functions  are  calculated  using  *-2 

CT2  =  e*2  h  -2/3  (z/h)"4/3  B  (7a) 

CQ2  =  M*2h-2/3  (Z/h)-4/3AB  (7b) 

CTQ  =  rTQ  M*©*h"2/3  (z/h)"4/3  A1/2b  (7c) 

As  has  been  pointed  out  by  Panofsky  13,  Eq .  7  is  an  extension  of  Eq .  3  for  |  Z/LJ  >>1. 
Thus,  the  matching  of  the  surface  layer- free  convection  layer  scaling  requires  that  B  =  2.7, 
which  has  been  verified  by  measurements  ** . 

Transition  Layer.   For  Z  >  0.5h,  Eq .  7  begins  to  breakdown.   The  actual  height  limit  of 
the  validity  of  free  convective  scaling  depends  on  conditions,  most  notably  the  value  of  h. 
Basically,  the  greater  h,  the  larger  the  value  of  Z/h  that  Eq .  7  can  be  used.   The  behavior 
of  C<p2  in  the  transition  layer  has  been  examined  using  an  additional  dimensionless 
function  1 4>  15  G(Z/h) 

CT2  =  e*2  h-2/3  (z/h)"4/33  G(Z/h)  (8) 

Three  examples  of  G(Z/h)  are  shown  in  Fig.  6,  where  it  is  clear  that  G(Z/h)  is  not  a 
universal  function.   Since  turbulence  in  the  transition  layer  is  influenced  by  the  transport 
of  turbulence  from  the  interfacial  layer,  it  is  likely  that  scaling  in  the  region  will 
require  "matching"  with  the  entrainment  scaling. 

Interfacial  Layer.   Turbulence  at  the  top  of  the  mixed  layer  erodes  the  interface 
separating  the  turbulent  boundary  layer  and  the  quiescent  ,  stable  air  of  the  free 
troposphere.   The  rate  of  this  process,  which  is  called  entrainment,  depends  upon  the  amount 
of  turbulence  available  and  the  amount  of  work  the  turbulence  must  do  to  erode  the  boundary. 
The  amount  of  turbulence  is  generally  a  function  of  the  surface  fluxes  and  the  mixing 
velocity,  W* .   The  barrier  the  turbulence  works  against  is  primarily  a  function  of  the 
bouyancy  jump,  A9V  and  the  lapse  rate,  TT,  of  the  upper  layer.   The  mixing  in  of  the 
warm,  dry  air  from  above  the  inversion  creates  considerable  fluctuation  in  temperature  and 
water  vapor  with  a  resultant  large  increase  in  0^2,  Cq^  ,  and  Ca- 
using Deardorff's  1*  model  of  the  interfacial  region,  Wyngaard  and  LeMone  *■*   have 
developed  a  model  to  predict  the  structure  functions 

CT2  =  A9V  (1  +  a)rThl/3  (We/W*)FTT  (9a) 

CQ2  =  4(AQ)2(l  +  ct)rfh1/3(We/W*)/A0v  (9b) 

CTQ  =  AQ(1  +  ajr-rh1/3  (We/W*)FTQ  (9c) 

where  AQ  is  the  jump  in  water  vapor  density  (remember  that  AQ  is  negative!),  We  is  the 
entrainment  rate  (m/s)  ,  F^q  and  F^t  are  functions  that  relate  virtual  temperature 
fluctuations  to  the  temperature  fluctuations  *-2m 

The  final  step  in  obtaining  a  structure  function  formulation  is  to  select  a 
parameterization  of  We .   Wyngaard  and  LeMone  chose  the  second  order  entrainment  model 
equilibrium  rate1^ 

We  =  0.8  Qov/((l  +  o)  TTh)  (10) 

and  obtain  excellent  agreement  with  NCAR  aircraft  data.   A  much  more  common  assumption  is 
that  17 

We  =  0.2QOV/A9V  (11) 
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In  fact,  the  major  uncertainty  in  the  application  of  Eq .  9  is  the  entrainment  energy 
(bouyancy,  surface  shear,  inversion  shear,  and  atmospheric  waves)  and  dozens  of  entrainment 
relationships^  exist.   Another  problem  is  the  appropriate  value  of  the  interface 
thickness,  a.   Typical  values  of  a  range  from  0.1  to  0.5  with  an  average  of  0.3.   Also, 
since  a  represents  an  ensemble  average  property,  a  single  measurement  (by  a  radiosonde,  for 
instance)  may  not  give  an  accurate  estimate.   Mote  that  if  Eq .  10  is  used,  a  and  T^, 
drop-out  of  the  structure  function  calculation  ( Eq .  9).   Another  source  of  error  in  the 
Wyngaard  and  LeMone  model  is  the  assumption  that  the  average  rate  of  dissipation  of 
turbulent  kinetic  energy,  £^,  in  the  interfacial  layer  is  simply  one  half  the  free 
convection  prediction 


£ 


-  1/2  £fc  =0.2  gQov/T  (12) 


Since  the  CN2  prediction  is  proportional  to  £"1/3,  an  underestimate  of  £j_  would 
lead  to  an  overestimate  of  CN2 ,   Thus,  wind  shear  at  the  inversion  not  only  leads  to  an 
underestimate  of  We  but  also  to  an  underestimate  of  £  j_  (note  that  these  two  errors  tend 
to  compensate).   A  summary  of  four  cloud-free  measurements  in  the  interfacial  layer  made 
during  the  MAGAT  experiment^  is  shown  in  Table  2. 


Table    2 .       Boundary    layer   and 

inter: 

Facial    laver 

c 

iata    from 

f< 

Dur 

MAGAT    aircraft    orofiles 

Date                                            4/3 

5/4 

5/4 

5/7 

Time,    PST                               1610 

102  4 

1201 

10  43 

h,    m                                           320.00 

300.00 

330.00 

200.00 

a                                                           .35 

.5 

.5 

.5 

A9V,K                                              6.5 

11.00 

9.00 

7.00 

AQ,g/m3                                        -4.4 

-5.2 

-5.2 

-2.00 

QOV/    Km/s                                       .027 

.022 

.02 

.02 

rT,    K/km                                      9.00 

10.00 

15.00 

9.00 

W*<    m/s                                            .65 
Bi1/3,    m2/3/s                            .13 

.6 
.06 

.7 

.12 

.5 
.067 

(0.2g   Qov/T) 1/3,m2/3/s       .056 

.053 

.052 

.052 

FT^                                                     .90 
CT2i      K2/m2/3                            3-  4   x 

.78 

.86 

.66 

10- 

-3 

7.3    x 

10" 

4 

3.3 

X 

10" 

-3 

6.0    x 

io-3 

Cq2i,     (g/m3)2/m2/3              9.6    x 

io- 

-3 

1.8    x 

10" 

3 

6.  4 

X 

10" 

-3 

1.9    x 

io-3 

Table  3  compares  the  measured  structure  functions  with  the  Wyngaard  and  LeMone  model 
using  Eq .  10  and  Eq .  11. 

Table  3.   Comparison  of  MAGAT  interfacial  structure  function  data 

with  entrainment  model  predictions 

Source  <CTZ>,  K27m2/3     <Cq2> ,  (g/mJ) 2/m2/3 

Measurements  3.5  x  j.q-3  4.9  x  ]_q-3 

We  =  0.3  Qov/(d  +  a)TTh)      4.5  x  in-3  5.6  x  IO*3 

We  =  0.2  Qov/ASv  0-6  x  IO-3  0.7  x  IO"3 

Free  troposphere 

Except  in  regions  of  clear  air  turbulence  (CAT),  Cm2  in  the  free  troposphere  is  quite 
small.   Typical  values  of  C-j.2  are  on  the  order  of  10-'  ^2/m2/3(  leading  to 
C*j2  %  5  x  10--9  m-2/3  at  Z  =  5km.   The  decrease  of  atmospheric  pr'essure  with 
increasing  altitude  is  another  factor  that  reduces  Cj,j2  relative  to  C-p2  (Fig.  7). 
There  are  measurements  in  this  region  20,21,22,23,2  4  jjUt  there  are  no  physical  models. 
There  are  several  phenomenological  models  21,25  that  relate  Op2  to  local  temperature 
vertical  gradients.   A  good  review  of  the  subject  is  given  by  Hall2^. 

Aerosol  modeling 

Given  the  validity  of  the  micrometeorological  relationship  of  molecular  and  aerosol 
extinction  to  actual  optical  extinction,  the  optical  model  is  considered  to  be  equivalent  to 
the  meteorological.   For  molecular  extinction,  the  pressure,  temperature  and  relative 
humidity  are  obtained  within  the  framework  of  the  MBL  structure  as  depicted  in  Fig.  1. 
Thus,  surface  layer  measurements  of  temperature,  relative  humidity,  plus  the  height,  h,  are 
sufficient  to  calculate  the  molecular  extinctions  (using  LOWTRAN)  throughout  the  MBL.   If, 
in  addition,  a  surface  layer  measurement  of  the  aerosol  spectrum  is  available,  then  it  is 
also  possible  to  calculate  the  aerosol  extinction  profile  throughout  the  MBL  using  an 
aerosol  mixed  layer  structure  which  we  shall  describe.   Although  temperature  and  humidity 
are  easily  available  over  the  ocean,  routine  measurements  of  aerosol  spectra  are  most 
certainly  far  in  the  future.   Thus,  it  is  necessary  to  develop  a  model  which  allows 
estimates  of  the  aerosol  spectrum  from  routinely  available  oceanic  meteorological  data.   The 
remainder  of  this  section  will  be  devoted  to  the  aerosol  problem  since  the  molecular  is 
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considered  to  be,  by  comparison,  more  than  adequately  solved. 

Aerosol  basics 

Before  attacking  the  aerosol  MBL  problem,  it  is  necessary  to  establish  some  general 
background  in  marine  aerosols. 

Relative  humidity.   It  is  well  known  that  aerosol  particles  over  the  ocean  are  usually  in 
the  form  of  dissolved  or  coated  droplets.   The  size  of  these  droplets  changes  as  liquid 
water  is  absorbed  or  lost  in  response  to  its  relative  humidity  variations2'.   For 
example,  let  rQ  be  the  particle  radius  at  some  reference  water  vapor  saturation 
SQ(S  =  RH/100)  and  r  be  the  particle  radius  at  ambient  saturation  S,  then 

r  =  r0  g(S)  (13) 

where  (assuming  SQ  =0.8) 

g(S)  =  0.81  exp(0.066  S/(1.058  -  S))  (14) 

This  relation  is  particularly  useful  because  one  can  eliminate  variations  of  the  aerosol 
spectrum  due  solely  to  humidity  variations  by  dealing  with  the  spectrum  as  it  would  appear 
at  the  reference  humidity.   Thus,  if  N(r)  is  the  ambient  aerosol  number  density,  we  can 
transform  N(r)  to  the  reference  saturation  spectrum  N(r0), 

N(r/g(S))  =  N(rQ)  =  N(r)  g(S)  (15) 

Since  some  of  the  results  to  be  described  are  in  aerosol  volume  format,  V(r)  =  4/3irr^  N(r), 
the  equivalent  volume  relation  should  be  considered 

V(r/g(S))  =  V(r0)  =V(r)/g2(S)  (16) 

Aerosol  spectrum  components.   Over  the  ocean,  the  aerosol  spectrum  at  small  sizes 
(r  <  1.0  )in)  is  dominated  by  particles  of  nonlocal  origin  (continental  background)  while  the 
larger  sizes  (r  >  1.0  ]xm)    are  primarily  sea  salt  droplets  produced  by  whitecaps.   Since 
particles  with  radius  less  than  0.1  ym  do  not  contribute  significantly  to  extinction  at 
visible  and  IR  wavelengths,  the  continental  aerosols  can  be  reasonably  represented  by  the 
Junge  distribution 

Vc(r0)  =  A/rQ  (17) 


If  we  define  Vg(r0)  as  the  sea  salt  component,  then  the  total  aerosol  spectrum  in  the 
MBL  is 

V(rQ)  =  Vs(r0)  +  Vc(r0)  .  (18) 

Accurate  broad  band  aerosol  spectrum  measurements  can  be  used  to  determine  the 
continental  influence  by  assuming  that  at  rQ  =  0.1  pin,  Vs(0.1)  =  0  and  therefore 

A  =  0.1  Vc(0.1)  =  0.1  V  (0.1)  (19) 

This  assumption  is  corroborated  by  aerosol  data  taken  in  the  North  Atlantic  (JASIN  exper- 
iment) that  has  been  averaged  in  ensembles  based  on  wind  speed  (Fig.  8).   The  lack  of  wind 
speed  dependence  of  the  aerosol  volume  at  rQ  =  0.1  pm  demonstrates  that  these  particles 
are  primarily  of  continental  origin.   The  average  value  of  A  obtained  was  1.0.   The  claim 
that  A  is  an  index  of  continental  influence  is  further  validated  by  a  comparison  of  NPS 
values  of  A  with  atmospheric  Radon  activity  obtained  by  NRL28  during  CEWCOM-78  (Fig.  9). 
Given  a  measurement  of  A,  one  can  now  calculate  the  sea  salt  aerosol  spectral  component 

Vs(r0)  =  V(r0)  -  A/r0  (20) 

One  note  of  caution:  in  order  to  be  consistent  with  the  literature,  we  will  use  aerosol 
coefficients  A  and  B;  these  bear  no  relation  to  A  and  B  in  the  section  on  turbulence. 

Aerosol  MBL  structure 

The  necessity  for  the  partitioning  of  aerosol  components  becomes  clear  when  one  considers 
the  different  particle  sources  they  represent.   Because  sea  salt  particles  are  produced  by 
whitecaps  at  the  sea  surface,  they  tend  to  remain  trapped  within  the  MBL.   Thus  we  assume 
Vg  =  0  for  Z  >  h.   Continental  aerosols  are  assumed  to  be  distributed  with  approximately 
constant  mixing  ratio  both  above  and  below  the  MBL.   The  vertical  structure  being  proposed 
is  depicted  in  Fig.  10.   Recall  that  we  are  describing  the  aerosol  density  at  reference 
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saturation,  S0  ■  O.R.   Since  the  saturation  tends  to  increase  with  increasing  altitude  in 
the  MBL  (this  is  consistent  with  Fig.  1  where  ©v  and  dewpoint  are  nearly  constant, 
therefore  T  must  decrease  roughly  one  degree  per  hundred  meters),  the  actual  measured 
spectrum  (or  extinction  coefficient)  will  tend  to  increase .  2<^<  30   Conversely,  if  one 
takes  a  measured  vertical  profile  of  aerosol  volume  and  removes  the  humidity  dependence 
using  Eg.  16,  then  the  reference  humidity  aerosol  volume  is  nearly  height  independent  in  the 
MBL  (Fig.  11).   Note  the  aerosol  volume  decreased  rapidly  above  the  MBL  indicating  that,  in 
this  case,  most  of  the  volume  in  the  MBL  is  due  to  sea  salt  particles. 

Surface  layer  gradient.   One  issue  not  covered  in  the  structure  shown  in  Fig.  10  is  the 
near  surface  vertical  gradient.   Since  continental  aerosols  are  not  surface  produced,  they 
do  not  have  significant  near  surface  gradients.   In  the  case  of  sea  salt  aerosols,  the 
surface  gradient  is  produced  by  a  near  balance  of  turbulent  upward  transport  of  particles 
and  the  gravitational  (Stokes  fallout)  downward  transport . 31 • 32   For  particles  in  the 
0.1  um  to  15  ym  range,  this  gradient  is  very  small  at  heights  greater  than  a  few  meters  and 
probably  not  worth  further  consideration  (except  for  long  wavelength,  near  surface  IR 
applications) . 

Scaling  law  perspective.   The  development  of  physical  scaling  equations  for  aerosols  is 
somewhat  simpler  than  for  C^2  because  one  is  not  modeling  the  fluctuations  but  rather 
the  mean  aerosol  profile.   Thus,  the  aerosol  model  is  roughly  analogous  to  the  mean  profile 
of  water  vapor  -  there  is  surface  production  (evaporation)  and  continental  background  (water 
vapor  advected  over  the  ocean)  which  is  present  above  and  below  the  MBL.   Unlike  the  C^-2 
case,  however,  the  surface  layer  aerosol  density  cannot  be  "a  priori"  calculated  without 
consideration  of  mixed  layer  processes  (this  is  also  true  of  the  mean  water  vapor  density). 
This  will  be  discussed  in  later  sections. 

The  Wells-Munn-Katz  (WMK)  model 

The  present  state-of-the-art  in  bulk  aerosol  models  is  represented  by  the  WMK 
model33,34  where  the  continental  component  is  the  Junge  type  (Eq.  17). 

Nc(r)  =  3/r04  (21a) 

B  =  3A/(4ir)  (21b) 

and  the  sea  salt  component  has  a  wind  speed  and  height  dependence  of  the  gamma  function 
form29 

Ng(r)  =  F(r0,u)  exp(-Z/hQ)  (22) 

(rQ  =  r/g(S)  as  in  Eq .  13  and  Eq .  14). 

The  exponential  decay  with  scale  height,  hQ  =  800m,  though  in  clear  disagreement  with 
the  well-mixed  assumption  (Fig.  11),  is  probably  adequate  for  some  applications  (satellite 
remote  sensing,  for  example).   Although  there  are  ongoing  efforts  to  improve  the  sea  salt 
modeling  function35r  j_t  j_s  certainly  worthwhile  to  examine  in  detail  the  near  surface 
performance  of  the  WMK  model. 

WMK  model  evaluation  using  JASIH  data.   This  evaluation  is  basically  a  comparison  of  WMK 
model  predictions  of  aerosol  extinction  based  on  measurements  of  relative  humidity  and  wind 
speed  with  extinctions  obtained  from  measurements  of  aerosol  spectra.   A  twenty  minute 
averaging  period  was  used  for  RH  and  aerosol  spectra.   The  wind  speed  used  in  Eq.  22  was  an 
average  over  the  12  hours  preceding  the  aerosol  data  in  order  to  take  into  account  the  long 
response  times  required  for  sea  salt  aerosol  generation.   This  comparison  (Fig.  12)  asks 
"given  a  wind  speed  and  relative  humidity,  how  well  does  the  model  predict  the  observed 
extinction"?   Clearly,  the  WMK  model  predicts  very  well  on  the  average  with  a  rather  large 
standard  deviation  (about  half  an  order  of  magnitude,  or  a  factor  of  3 ) .   A  considerably 
more  stringent  comparison  asks  "given  an  observed  extinction,  how  well  does  the  model 
predict  this  extinction"?   This  comparison  (Fig.  13)  is  considerably  less  favorable, 
particularly  in  the  extremes  of  good  visibility  and  low  visibility  conditions.   A  certain 
amount  of  disagreement  of  this  type  is  expected  when  comparing  quantities  that  are  subject 
to  experimental  error.   However,  we  found  that  a  considerable  part  of  the  discrepancy  is  due 
to  overestimation  of  the  continental  aerosol  component  (WMK  assume  3  =  1.7  while  JASIN  gives 
B  =  0.2  4).   If  we  change  3  to  0.2  4  and  repeat  the  comparison  (X's  in  Fig.  13),  we  get  some 
improvement.   Furthermore,  there  is  considerable  correlation  between  visible  extinction  and 
the  amount  of  continental  aerosol;  the  high  visibility  conditions  have,  relatively,  much 
smaller  continental  aerosol  densities  while  the  low  visibility- conditions  have  relatively 
greater  continental  aerosol  densities.   The  disagreement  at  high  extinction  values  is  due  to 
the  inaccuracy  of  the  relative  humidity  measurement  which  is  critical  under  near  fog  and 
heavy  haze  conditions,  where  RH  =  100%. 
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It  is  important  to  realize  that  the  WMK  model  is  an  average  continental,  equilibrium  sur- 
face generation  model.   Variations  of  the  continental  aerosol  were  more  than  an  order  of 
magnitude  during  JASIN.   The  sea  salt  component  is  also  subject  to  considerable  deviations 
from  equilibrium  (this  is  one  of  the  reasons  for  the  12  hour  average  wind  speeds).   The  dev- 
iations from  average  and  equilibrium  values  are  more  important  for  operational  usage  and 
somewhat  less  important  for  climatological  and  spectrum  evaluation  usage.   An  example  where 
both  usages  are  affected  is  the  estimation  of  IR  extinction  using  visible  extinctions  (visi- 
bility observations)  and  empirical  assumptions  about  their  relationship.   We  have  plotted 
the  visible  to  IR  extinction  ratio  in  Fig.  14  with  lines  indicating  the  WMK  model  predic- 
tions using  two  different  continental  aerosol  coefficients.   Since  the  continental  aerosol 
coefficients  are  correlated  to  visibility,  the  correct  ratio  depends  not  only  on  wind  speed 
but  also  the  visibility  observation. 

Although  we  have  dealt  with  the  stochastic  properties  of  ensemble  averages  of  the  aero- 
sols, variations  about  the  average  are  not  necessarily  random  but  are  primarily  due  to 
changes  in  synoptic  and  mesoscale  weather  patterns.   In  the  case  of  the  continental  aerosol 
component,  this  is  basically  a  question  of  air-mass  history.   In  the  case  of  the  sea  salt 
aerosol,  it  is  a  question  of  changing  surface  generation  rates  (wind  speed)  and  the  produc- 
tion, removal  and  mixing  mechanisms  in  the  marine  atmospheric  boundary  layer.   Since  the 
surface  generated  aerosols  are  quickly  mixed  verically  to  fill  the  boundary  layer  up  to  the 
capping  inversion,  rapid  changes  in  the  inversion  height,  h,  will  be  reflected  in  changes  in 
the  sea  salt  aerosol  density  and,  therefore,  the  extinction  coefficient.   In  Fig.  15  we  can 
see  that  the  fractional  variations  in  h  are  highly  correlated  with  variations  in  the  visible 
extinction.   The  correlation  with  10.6  um  extinction  is  considerably  less  because  the  large 
size  aerosols  (which  are  heavier  contributors  to  IR  extinction)  reach  equilibrium  more 
quickly  after  changes  in  surface  conditions. 

The  aerosol  mixed-layer  model 

The  sea  surface  is  a  continuous  source  of  sea  salt  aerosols  in  the  marine  boundary  layer. 
These  surface  produced  aerosols  can  be  characterized  by  a  surface  flux  spectrum,  Fs(rQ), 
which  represents  the  volume  of  aerosol  per  particle  radius  interval  produced  per  square 
centimeter  of  ocean  surface  each  second  as  a  function  of  aerosol  particle  radius.   This 
quantity  is  a  function  of  wind  speed. 

The  continuous  production  of  sea  salt  aerosols  is  balanced  by  several  removal  mechanisms. 
One  obvious  mechanism  is  the  loss  of  particles  as  they  fall  back  to  the  surface.   This  set- 
tling under  gravity  is  called  "Stokes  fallout"  and  is  characterized  by  the  Stokes  velocity. 
The  particles  are  transported  vertically  by  turblence  in  the  marine  boundary  layer  and 
maintained  at  a  uniform  density  throughout  the  mixed  layer.   The  growth  of  the  height  of  the 
layer  constitutes  another  loss  mechanism  called  "entrainment " ,  which  was  discussed  in  the 
interfacial  turbulence  section  of  this  paper.   The  final  loss  mechanism  is  "rainout"  which 
occurs  when  the  particles  become  condensation  nuclei  in  the  formation  of  clouds. 


Given  the  surface  flux  spectrum  and  a  parameterization  of  the  removal  mechanisms,  one 
could  predict  evolutions  of  the  aerosol  density  spectrum  by  the  following  equation-^ 


hdVs/dt  =  Fs  -  (We  +  WS)VS  (23) 

where  Wg  is  the  Stokes  velocity 

Ws  =  1.57  x  lO-2  (1  +  1.2  (0.9l/g(S))3)ro2g2(S)  (24) 

where  rQ  is  in  um,  Ws  is  in  cm/sec.   In  Eq.  23  we  have  left  out  the  cloud  formation 
removal  mechanism  because  it  acts  on  a  much  longer  time  scale.   One  assumption  implicit  in 
Eq.  2  3  is  that  the  convective  mixing  velocity,  W* ,  is  much  larger  than  Ws .   Since  W* 
is  typically  1  m/sec  while  Wg  is  on  the  order  of  0.1  m/sec  for  the  largest  particles 
(r  =  15  um)  considered  here,  the  assumption  is  reasonable. 

The  surface  flux  spectrum  Fs(r0),  which  is  a  major  factor  in  Eq.  23,  is  not  well 
known.   Some  preliminary  estimates  are  available  from  laboratory  measurements^.   In  an 
effort  to  obtain  improved  information  about  Fs(r0),  the  CEWCOM-7  8  data  have  been  used  to 
evaluate  the  surface  flux. ^°      Since  the  other  terms  of  the  expression  were  measured,  one 
simply  calculates  Fs  using  a  rearrangement  of  the  terms.   Thus,  the  final  expression  for 
the  flux  calculation  is, 

Fg   =   hdVg/dt   +    (Wg   +   Wg)Vg  (25) 

A  period  was  chosen  from  CEWCOM-78  where  all  data  were  available,  the  synoptic  conditions 
were  fairly  constant,  the  wind  speed  was  reaonably  constant  and  a  good  mixed  layer  was 
present  (5/20-5/21).   Because  of  the  reasonably  constant  wind  speeds  during  the  20  hour 
CEWCOM-78  period,  we  were  able  to  improve  the  statistical  certainty  for  the  F£  calculation 
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by  combining  all  the  data  from  the  period,  assuming  that  the  flux  would  be  a  reasonable 
representation  for  U  =  9  m/sec  wind  soeed.   Thus,  we  now  have  the  surface  flux  spectrum  at  a 
single  wind  speed.   In  order  to  estimate  the  flux  at  other  wind  speeds,  we  note  that  the 
right  hand  side  of  Eq .  25  is  nearly  independent  of  wind  speed  for  equilibrium  conditions. 
Therefore,  the  flux  at  one  wind  speed  can  be  related  to  the  flux  at  other  wind  speeds  if  the 
equilibrium  volume  spectra  are  known: 

F3(Ul)    =  FS(U2)  VS(U1)/VS(U2)  (26) 

We  have  available  from  JASIN  a  large  set  of  ensemble  averages  of  aerosol  volume  spectra 
at  different  wind  speeds  (Fig.  16).   It  is  a  simple  matter  to  apply  this  data  to  Eq.  26, 
using  the  CEWCOM-78  aerosol  flux  and  equilibrium  spectrum  to  generate  the  surface  volume 
flux  spectra  as  a  function  of  wind  speed  (Fig.  17). 

Equilibrium.   During  periods  of  fairly  constant  boundary  layer  conditions,  the  aerosol 
spectrum  may  be  in  a  state  of  equilibrium  (that  is,  dVs(rQ)/dt  =  0).   Under  these 
conditions,  the  average  equilibrium  spectrum  is 

V3(r0)  =  F3(r0)/(We  +  Ws)  (27) 

— 

^s^ro^'  which  is  a  function  of  wind  speed,  humidity  and  entrainment,  is  analogous  to  the 
WMK  sea  salt  spectrum  which  is  also  for  an  "average"  condition.   Deviations  of  Vg(r0) 
from  Vs(r0)  are  described  by  Eq .  23. 

Non-equilibrium.   For  time  periods  of  a  few  hours,  the  aerosol  spectrum  may  not  be  in  a 
state  of  dynamic  equilibrium.   If  we  rewrite  Eq.  23  in  the  following  form: 

dVg/dt  +  Vs/Tp  =  Fs/h  (28) 

where  the  time  constant,  Tp»  is 

Tp  =  h/(We  +  Ws)  (29) 

then  we  see  the  analogy  of  aerosols  and  a  capacitor  charged  by  an  applied  "voltage",  Fs , 
through  a  "resistance",  (We  +  Ws)~l.   In  this  analogy,  the  "capacitance",  is  h. 

The  response  time  of  the  aerosol  density  is  a  strong  function  of  particle  radius  because 
Wg  <v  r2.   Values  of  tp  for  S  =  0.8  and  h  =  400  m  are  given  in  Table  4. 

TABLE  4   The  equilibrium  time  constant  for  aerosols 
at  S  =  0.8,  h  =  400  and  We  =  0. 4  cm/s 

r.  um   0.5      1  5  10      15 


Tp,  hours      28       22     11       3.5     0.5 


The  boundary  layer  mixing  time,  T^,  which  represents  the  time  required  for  changes  in 
aerosol  density  to  be  evenly  distributed  throughout  the  marine  layer  is: 

Tm  =  h/W*  {30) 

For  the  CEWCOM-78  analysis  period  (h  =  400  m  and  W*  =  0.6  m/s)  we  find  Tm  =  .16  hours. 
Therefore,  short  term  variations  in  mixing  volume  (dh/dt)  will  lead  to  changes  in  the 
aerosol  density  because  the  production  response  time  is  much  slower  than  the  mixing  time. 
Thus,  for  time  periods  on  the  order  of  one  hour,  changes  in  the  aerosol  density  (dVs/dt) 
will  be  highly  correlated  with  the  mixing  volume  term  (h_i  Vs  dh/dt) .   The  effect  will 
be  particularly  noticeable  for  smaller  particles  and  is  the  explanation  for  Fig.  15. 

Discussion 

Given  the  great  difficulty  of  direct  optical  measurement  of  the  propagation  properties 
(CN2  an(j  extinction  vertical  profiles)  of  the  MBL,  it  is  necessary  to  estimate  these 
properties  from  more  practical,  routine  atmospheric  data.   This  is  accomplished  by  relating 
the  optical  properties  to  micrometeorological  quantities  (turbulence,  molecular  constituents 
and  aerosol  spectra)  where  the  atmospheric  physical  foundation  exists  for  the  development  of 
reasonably  simple  scaling  models.   The  first  step  in  this  process  was  to  verify  the  opti- 
cal-meteorological relationships  in  the  marine  regime.   Although  Cnj2  an<^  extinction  were 
considered  separately,  both  models  were  cast  in  the  same  framework  of  a  well  mixed  marine 
layer  capped  by  an  entraining  inversion;  the  structure  of  the  MBL  i3  a  key  element  in  model- 
ing of  the  height  dependences  above  the  near  surface  region.   Since  the  mixing  and  entrain- 
ment process  (also,  the  turbulence)  are  driven  by  the  surface  fluxes,  the  key  element  of 
this  effort  was  the  development  and  verification  of  the  bulk  aerodynamic  method  which  allows 
calculation  of  the  fluxes  and  scaling  parameters  from  routine  surface  observations  (sea 

B-9 


surface  temperature,  air  temperature,  relative  humidity  and  wind  speed).   A  measurement  for 
estimation  of  the  mixed  layer  height,  h,  will  permit  calculation  of  the  profiles  for  most  of 
the  remaining  (non-surface)  MBL.   Both  CN2  (near  the  inversion)  and  aerosol  spectral 
structure  models  require  knowledge  of  the  entrainment  rate  which  necessitate  at  least  a 
radiosonde  estimate  of  the  temperature  and  humidity  discontinuities  at  the  inversion  and  the 
free  troposphere  temperature  lapse  rate.   The  aerosol  model  also  requires  a  specification  of 
the  continental  component.   At  present  this  is  an  assumed  climatological  mean  but  it  could 
be  related  to  visibility  observations  or  air-mass  trajectory. 

The  accuracy  of  these  models  is  very  difficult  to  access.   In  the  case  of  turbulence,  the 
accuracy  of  the  input  data  is  critical.   For  aerosols,  the  variations  of  continental  lerosols 
and  open  ocean  mixed  layer  effects  are  simply  not  very  well  known.   Based  solely  on  our 
experience,  a  best  guess  of  the  RMS  error  of  the  turbulence  model  is  about  a  factor  of  two 
(average  over  the  MBL);  for  extinction  the  error  is  about  a  factor  of  three. 
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Figure  1   Typical  structure  of  the 
marine  atmospheric  boundary  layer  in 
terms  of  water  vapor  density,  0,  and 
virtual  potential  temperature,  9V. 


B-ll 


C*   Turbulence 

Figure  2   Refractive-index  structure 
function,  Cnj2,  determined  from 
direct  optical  measurement  versus  C^2 
determined  from  -licrothermal  point 
measurements  of  C<r<2  (Eq.  1)  assuming 
Cq2   and  CTq  are  negligible.   The 
dashed  lines  represent  factors  of  2 
disagreement. 
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Figure  3   Aerosol  extinction 
coefficient,  a ,  determined  from  direct 
optical  measurement  versus  a  determined 
from  point  measurements  o^  aerosol 
spectral  density,  N(r) ,  and  Eq .  2. 


IV!  AG  AT  - 

80 

- 

1.06  m 

- 

0        / 

a      o   _/ 

"V. 

9 

*>/ 

\s    o      °°    ° 

9S; 

r     9 

ffjr 

o     o 

- 

o'J, 

o 

3   £&£ 

0 

Jfao 

l 

QyOO 

<*?/* 

I 

I 

i          i          i          i 

100 

1 

10 

_ 

43 

f 

T.2 

i 
313 
S 

/ 

1 

-4 

12 

29 

/84 

1 

| 

i                i 

i 

I 

6 

-4 

-2 

D                           2 

e 


Figure  4a   Dimensionless  temperature 
structure  function:   overwater 
measurements  (points)  and  overland 
empirical  function7  f(Z/L)  (line). 


Figure  4b   Dimensionless  water  vapor 
structure  function:   overwater 
measurements  .(points)  and  empirical 
function7  A  f(Z/L)  (solid  line). 


B-T2 


Ch2 

Opipcal 


kT 


CN3  Bulk 


Figure  5   Similar  to  Figure  2  but  the 
horizontal  axis  is  C^  determined 
from  bulk  model  calculations  (Eg.  1  and 
Eq.  3). 
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Figure   7 
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Figure  6   Average  mixed  layer 
dimensionless  temperature  structure 
function.   The  solid  line  is  Frisch  and 
Ochs14  (h=300m) .   The  dashed  line  is 
Kaimal  et  al10  (h=1700m)  and  the 
points  are  from  Fairall  et  al5-3 
(h=»700m) . 
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Figure  8   Ensemble  average  total  aerosol 
volume  spectra  from  JASIM.   The  number 
to  the  right  of  the  spectrum  is  the  wind 
speed  category  in  m/sec. 
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Figure  9   Time  series  of  atmospheric 
222fcn  activity  (solid  line,  Larsen, 
Kasemir  and  Bressan28)  an(j 
continental  .aerosol  coefficient  (dashed 
line) ,  A. 
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Figure    10      Schematic    representation   of 
the   height   dependence  of  continental 
(Vc)    and    sea-salt    (Vs)    aerosols    at 
S=S«. 
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Figure  11   MAGAT  aircraft  measurements 
of  the  vertical  profiles  of  relative 
humidity  (RH),  ambient  aerosol  total 
volume  (V)  and  reference  (3=S0) 
aerosol  total  volume  (V^).   The 
aerosol  volume  is  in  um^/cm  . 
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Figure  12   Aerosol  extinction 
coefficient  (circles)  and  WMK  model 
predictions  (line)  at  X   =  10.6  ym  and 
RH  =  87%  as  a  function  of  12  hour 
average  wind  speed. 


10" 


h 

II 

x 

LU 


10 


icr3Li 


A»10.6um 

o  B=1.7  12  hr  u  avg 
x  B=0.24  12  hr  u  avg 
a  B  from  a  (0.49) 


JASIN 


x 

_      A 


10"3  10  '  10 

Measured  Extinction  (km"1) 

Figure    13      WMK   model    predictions    of 
aerosol    extinction   at    \   -    10.6    ym   as    a 
function   of   observed    extinction    (circles 
are    B   =    1.7,    X's    are    B    =    0.24  and 
triangle    are    3   obtained    from   average 
correlation  with   visibility) . 
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Figure  1 4   Ratio  of  aerosol  extinction 
coefficient  at  X  =  0.49  ym  to  that 
obtained  at  X   =  3.75  ym  and  10.6  ym.  The 
circles  are  the  JASIN  data  at  X    =  10.6  ym., 
the  triangles  are  the  JASIN  data  at 
X   =  3.75  ym.   The  dashed  line  is  WMK 
model  with  B  =  0.24  while  the  dotted 
line  is  3  =  1.7. 
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Figure  15   Fractional  changes  of  aerosol 
extinction  coefficient  versus  fractional 
changes  in  mixed  lay«=>r  height  h  for 
successive  4  hour  periods. 
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Figure  17  Ensemble  average  surface  flux 
spectra  deduced  from  Eq .  26.  The  number 
to  the  right  of  the  spectrum  is  the  wind 
speed  in  m/sec. 
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Figure  16   Ensemble  average  equilibrium 
sea  salt  volume  spectra  from  JASIN.   The 
number  to  the  right  of  the  spectrum  is 
the  wind  speed  category  in  m/sec. 
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